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With  reference  to  Part  II,  Section  V  (Pg.  107 — 144)  of  the  report,  the  following 
background  information  is  useful  to  the  reader  before  noting  the  corrections  presented 
here. 

Because  the  five-day  sample  of  constant  pressure  prognoses  and  analyses  supplied 
to  us  by  3rd  Weather  Wing  (3rd  WW)  was  in  a  grid  array  on  magnetic  tape,  the  grid  arrays 
had  to  be  read  and  printed  at  United  Aircraft  Research  Laboratories  (UARL)  using  a 
procedure  for  3rd  WW  grids  that  has  frequently  been  used  in  the  past.  For  some  un- 

•s 

determined  reason,  these  values  on  tape  were  read  and  printed  by  the  computer  in  an 
order  directly  opposite  to  what  the  real  order  was  supposed  to  be.  This  was  true  for  the 
analyses,  the  prognoses,  and  the  error  fields.  Thus,  it  has  no  effect  whatsoever  on  the 
evaluation  of  the  diagnostic  LMW  prediction  technique  described;  this  fact  cannot  be 
overemphasized.  The  only  effect  it  does  have  on  the  results  presented  in  Section  V  is 
that  the  hemispheric  geography  is  flipped  relative  to  the  analysis,  prognostic,  and 
error  isopleths  (e.g.,  North  America  should  be  where  central  Asia  is  and  vice-versa). 

The  prognoses,  analyses,  and  error  fields  do  not  change. 

The  incorrect  reading  and  printing  of  the  grids  (Figures  5  through  12,  Part  II, 

Section  V)went  unnoticed  because  raw  data  had  not  been  available  from  3rd  WW  to 
check  observed  station  values  against  the  analyses  (due  to  security  classification  in 
force  at  the  time  for  these  data  tapes).  Recent  acquisition  of  the  station  radiosonde 
observations  resulted  in  the  discovery  of  the  incorrect  orientation  of  the  grids. 

Because  of  this  incorrect  orientation,  the  following  changes  are  necessary  to  the 
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P.  108  §  8,  112, 
lines  3  and  4 


but  over  central  and  eastern  but  over  North  America 

Asia  and  the  western  Pacific  and  the  western  Atlantic 


P.  109,  second  full  H,  forecasts  occurred  over  central 
line  4  and  eastern  Asia 


forecasts  occurred  over 
North  America 


P.  109,  second  full  H,  and  over  an  area  near  Hawaii 
line  5 


and  over  north-eastern 
Africa 


P.  109,  second  full  H,  Over  the  North  Atlantic  Ocean 
line  7 


Over  eastern  Asia  and 
the  western  Pacific 


P.  111-144, Figures  The  base  map  should  be  rotated  180°  (e.g.,  what  is  now 
5—12  100°E  longitude  should  be  80°W  longitude). 


TECHNIQUES  FOR  OBJECTIVE 
HEMISPHERIC  ANALYSIS  AND  PREDICTION 
OF  THE  JET  STREAM 


David  B.  Spiegler 
John  T.  Ball 
et  al. 


September  1965 


433L  SYSTEM  PROGRAM  OFFICE 
ELECTRONICS  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
L.  G.  Hanscom  Field,  Bedford,  Mass. 


FOREWORD 


System  433L;  project  2.0;  task  2.3.  This  TR  has  been  prepared  for  United 
Aircraft  Corporation,  East  Hartford,  Conn.,  under  Subcontract  No.  15107  to 
Contract  No.  AF19(628)-3437,  by  The  Travelers  Research  Center,  Inc.,  250  Constitution 
Plaza,  Hartford,  Conn.  The  Research  Center’s  publication  number  is  7463-176. 

Robert  L.  Houghten,  Lt.  Colonel,  USAF,  is  Acting  System  Program  Director.  This 
report  covers  the  period  1  October  1964—31  July  1965,  and  was  submitted  for  approval 
on  7  September,  1965. 

The  authors’  names  not  listed  on  the  front  cover  are  Bernard  J.  Erickson  and 
Robert  J.  Ball. 

The  programming  support  for  the  analysis  and  prediction  techniques  described 
in  this  report  was  outstanding.  Mr.  Charles  Pike  of  UAC  had  prime  responsibility  in 
this  area  and  was  ably  assisted  by  Mrs.  Lydia  Rufleth.  The  authors  wish  to  thank 
them  for  a  job  that  was  extremely  well  done. 

Thanks  are  also  due  to  the  many  people  in  the  Meteorological  Analysis  and 
Prediction  Division  of  TRC  who  participated  in  various  phases  of  the  work— particularly 
to  Mr.  Keith  W.  Veigas  for  his  keen  interest  and  helpful  suggestions  during  discussions 
of  the  jet-stream  analysis  and  prediction  problems. 

The  National  Meteorological  Center  supplied  the  data  for  use  in  this  study. 


ABSTRACT 


This  report  comprises  two  parts:  Part  I  documents  an  objective  analysis  tech- 
nique  for  level  of  maximum  wind  (LMW)  parameters,  and  Part  II  presents  two  LMW 
prediction  techniques. 

The  objective  hemispheric  analysis  technique  incorporates  information  from 
data  areas  into  LMW  analyses  for  no-data  areas.  This  is  done  in  a  quantitatiye, 
objective,  and  consistent  manner  by  using  LMW  regression  equations  that  specify  the 
initial-guess  fields  for  the  analysis.  The  analysis  technique  locates  jet  cores  between 
grid  points,  and  generates  observations  along  these  cores  by  using  horizontal  jet  pro¬ 
files.  This  technique  produces  high-quality  analyses  that  compare  favorably  with 
subjective  analyses,  and  that  are  generally  consistent  with  one  another. 

The  first  prediction  technique  is  based  on  a  procedure  wherein  potential  pre¬ 
dictors  are  obtained  from  a  network  grid  oriented  with  respect  to  the  LMW  wind 
direction  at  the  predictand  point.  Prediction  equations  are  obtained  by  applying  the 
screening  regression  technique  to  predictors  selected  from  initial-state  LMW  param¬ 
eters  (and  their  12-hr  change)  and  wind  speed  values  at  constant  pressure  surfaces. 
Tests  on  independent  data  showed  that  these  equations  produce  prognoses  that  are  better 
than  persistence  for  all  polar  and  subtropical  jet-stream  categories. 

The  second  prediction  technique  applies  diagnostic  LMW  regression  equations 
(originally  derived  to  produce  the  initial-guess  fields  in  the  LMW  analysis)  to  constant- 
pressure -surface  numerical  prognoses.  A  five-day  data  sample  of  numerical  prognoses 
was  used  to  test  the  technique.  Evaluation  of  the  resulting  LMW  prognoses  indicated 
that  the  12-  and  24-hr  prognoses  were  generally  good  representations  of  the  observed 
LMW  fields,  and  that  the  36-  and  48-hr  prognoses  were  occasionally  good  over  some 
areas  of  the  hemisphere. 


REVIEW  AND  APPROVAL 

Publication  of  this  technical  report  does  not  constitute  Air  Force  approval  of  the 
report’s  findings  or  conclusions.  It  is  published  only  for  the  exchange  and  stimulation 
of  ideas. 


Acting  System  Program  Director 


iii 


TABLE  OF  CONTENTS 


PART  I 

Section  Title  Page 

I  INTRODUCTION  1 

II  LMW  SCREENING  REGRESSION  EXPERIMENTS  2 

III  METHOD  OF  ANALYSIS  AND  VERIFICATIONS  OF  LMW  PARAMETERS  8 

1.  Wind  direction  at  the  LMW  [W^(L)]  8 

2.  Wind  speed  at  the  LMW  [W  (L)]  8 

s 

3.  Verification  of  W  (L)  9 

s 

4.  Analysis  and  Verification  of  Z(L)  11 

5.  Analysis  and  Verification  of  S^^  and  S^  11 

IV  RESULTS  12 

6.  Results  of  LMW  Regression  Experiments  12 

7.  Verification  and  Evaluation  of  LMW  Analyses  22 

8.  Areal-mean -error-method  Statistics  22 

9.  Objective  (computer)  versus  Subjective  (man)  LMW  Analysis 

Comparison  32 

10.  Consistency  Between  LMW  Parameters  41 

V  CONCLUSIONS  AND  RECOMMENDATIONS  ,  52 

APPENDIX.  EQUATIONS  FOR  LMW-PARAMETER  FIRST-GUESS  FIELDS  57 

REFERENCES  60 


LIST  OF  ILLUSTRATIONS 

Figure  Title  Page 

1  Average  horizontal  profile  of  wind  speeds  in  the  jet  stream 

(taken  from  [2])  10 

2  Initial  guess  rms  errors  for  LMW  wind  speed  at  analysis  and  with¬ 
held  stations  for  12  observation  times  (area  where  W  >  100  knots)  25 

3(a)  Objective  LMW  wind  speed  analysis,  1200Z,  January  8,  1963  (knots)  26 


V 


Figure  Title  Page 

3(b)  Objective  LMW  wind  speed  analysis,  1200Z,  January  9,  1963  (knots)  27 

2 

4(a)  300-mb  analysis,  1200Z,  January  8,  1963  (10  ft)  28 

2 

4(b)  200-mb  analysis,  1200Z,  January  8,  1963  (10  ft)  29 

2 

5(a)  300-mb  analysis,  1200Z,  January  9,  1963  (10  ft)  30 

2 

5(b)  200-mb  analysis,  1200Z,  January  9,  1963  (10  ft)  31 

6  Rms  errors  of  wind  speed  at  the  LMW  for  three  wind-speed  categories 

for  12  observation  times  33 

7  Rms  errors  of  height  of  the  LMW  for  three  wind-speed  categories 

for  12  observation  times  34 

8  Rms  errors  of  shears  below  (a.)  and  above  (b)  the  LMW  for  three  wind- 

speed  categories  for  12  observation  times  35 

9(a)  Objective  and  subjective  jet  stream  analysis  comparison,  1200Z, 

January  1,  1963  (knots)  36 

9(b)  Objective  and  subjective  jet  stream  analysis  comparison,  1200Z, 

January  5,  1963  (knots)  37 

10(a)  Objective  and  subjective  height  of  LMW  analysis  comparison,  1200Z, 

January  1  1963  (10^  ft)  39 

10(b)  Objective  and  subjective  height  of  LMW  analysis  comparison,  1200Z, 

January  5,  1963  (10^  ft)  40 

11(a)  Objective  and  subjective  shear  below  the  LMW  analysis  comparison, 

1200Z,  January  1,  1963  (knot  10"^  ft)  42 

11(b)  Objective  and  subjective  shear  below  the  LMW  analysis  comparison, 

1200Z,  January  5,  1963  (knot  10“^  ft)  43 

12(a)  Objective  and  subjective  shear  above  the  LMW  analysis  comparison, 

1200Z,  January  1,  1963  (knot  10“^  ft)  44 

12(b)  Objective  and  subjective  shear  above  the  LMW  analysis  comparison, 

1200Z,  January  5,  1963  (knot  10"^  ft)  45 

13(a)  Objective  LMW  wind  speed,  initial  guess,  1200Z,  January  3,  1963 

(knots)  46 

13(b)  Objective  LMW  wind  speed,  final  pass,  1200Z,  January  3,  1963  (knots)  47 

13(c)  Objective  height  of  the  LMW  analysis,  1200Z,  January  3,  1963  (10^  ft)  48 

13(d)  Objective  shear  below  the  LMW  analysis,  1200Z,  January  3,  1963 

(knot  10”^  ft)  49 

13(e)  Objective  shear  above  the  LMW  analysis,  1200Z,  January  3,  1963 

(knot  10’3  ft)  50 


vi 


LIST  OF  TABLES 


Table  Title  Page 

I  Jet  stream  categories  and  their  definitions  4 

II  Relation  of  jet  category  definitions  to  height  of  the  LMW  6 

III  Selection  of  withheld  stations  on  the  basis  of  station  density  and 

LMW  wind  speed  11 

IV  Initial-guess  equations  for  LMW  parameters— dependent  and 

independent  data  results  for  seven  jet  stream  categories  13 

V  Percent  reduction  of  variance  associated  with  derived  and  non- 

derived  predictors  15 

VI  Overall  rms  errors  for  LMW  parameters  for  12  observation  times  23 


Vll 


TABLE  OF  CONTENTS 


PART  II 

Section  Title  Page 

I  INTRODUCTION  69 

II  DATA  PROCESSING  FOR  THE  PREDICTOR  NETWORK  TECH-  70 

NIQUE 

HI  THE  LMW  PREDICTOR  NETWORK  TECHNIQUE  72 

1.  Categorization  of  Predictand  Points  72 

2.  Predictor  Network  Grid  74 

3.  Predictands  75 

4.  Predictors  75 

IV  RESULTS  USING  THE  PREDICTOR  NETWORK  TECHNIQUE  79 

5.  Polar  Jet  Stream— Dependent  Data  Results  82 

6.  Subtropical  Jet  Stream— Dependent  Data  Results  86 

7.  Independent-data  Testing  86 

V  RESULTS  USING  THE  MODELING  PREDICTION  TECHNIQUE  107 

8.  Wind  speed  at  the  LMW  108 

9.  Height  of  the  LMW  109 

10.  Shear  below  the  LMW  110 

11.  Shear  Above  the  LMW  110 

VI  CONCLUSIONS  AND  RECOMMENDATIONS  145 

APPENDIX.  EQUATIONS  FOR  POLAR  AND  SUBTROPICAL  JET¬ 
STREAM  CATEGORIES  147 

REFERENCES  157 


viii 


LIST  OF  ILLUSTRATIONS 

Figure 

Title 

Page 

1 

Area  used  for  North  American  jet-stream  prediction  study 

71 

2 

Procedure  for  categorizing  NWP  grid  points  from  the  initial- 
state  position  of  the  jet  stream  core 

73 

3 

Predictor  network  grid 

75 

4  (a) 

W  (L)  analysis,  OOZ  December  17,  1963 

102 

(■b) 

\V  (Li  analysis,  i2Z  December  17,  1963 

8  * 

103 

(c) 

W  (Li  12-hr  prognosis  valid  12Z  December  17,  1963 
s 

103 

(d) 

W  (Li  analysis  OOZ  December  18,  1963 

S’- 

104 

(e) 

W  (Li  24-hr  prognosis  valid  OOZ  December  18,  1963 
s 

104 

(f) 

S^  24-hr  prognosis  valid  OOZ  December  18,  1963 

105 

(g) 

S  24-hr  prognosis  valid  OOZ  December  18,  1963 
a 

105 

(h) 

Z(L)  24-hr  prognosis  valid  OOZ  December  18,  1963 

106 

5  (a) 

W  (L)  analysis  OOZ  December  13,  1964 

S  ' 

111 

(b) 

W  (L)  12-hr  prognosis  valid  OOZ  December  13,  1964 
s 

112 

(c) 

W  (L)  24-hr  prognosis  valid  OOZ  December  13,  1964 
s 

113 

(d) 

W^(L)  36-hr  prognosis  valid  OOZ  December  13,  1964 

114 

(e) 

W^(Li  48-hr  prognosis  valid  OOZ  December  13,  1964 

115 

6  (a) 

W  (Li  analysis.  OOZ  December  14,  1964 
s  ' 

116 

(b) 

W  (Li  12-hr  prognosis  valid  OOZ  December  14,  1964 
s 

117 

^c) 

W  (L)  24-hr  prognosis  valid  OOZ  December  14,  1964 

118 

(d) 

W  (Lj  36-hr  prognosis  valid  OOZ  December  14,  1964 
s 

119 

(e) 

W^(L)  48-hr  prognosis  valid  OOZ  December  14,  1964 

120 

7  (a) 

ZfL)  12-hr  prognosis  (va’id  OOZ  December  13,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 
>  100  knots)]  ^ 

121 

(b) 

Z(L)  24-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 
>  100  knots)]  ® 

122 

(c) 

Z(Li  36-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 
>  100  knots'll  ^ 

123 

ix 

Page 


Figure 


Title 


7  (d)  Z(L)  48-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 

>  100  knots)] 

8  (a)  Z(L)  12-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 

>  100  knots)] 

(b)  Z(L)  24-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 

>  100  knots)]  ^ 

(c)  Z(L)  36-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 

>  100  knots)]  ^ 

(d)  Z(L)  48-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 

>  100  knots)]  ^ 

9  (a)  S,  12-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 

[(lorecast  minus  observed)  over  areas  where  W  (L)  observed  as 

>  100  knots)]  ^ 

(b)  Sia  24-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 

>  100  knots)]  ^ 

(c)  Sjj  36-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 
[(forecasting  minus  observed)  over  areas  where  W  (L)  observed  as 

>  100  knots)]  ^ 

(d)  S]3  48-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 

>  100  knots)]  ^ 

10  (a)  Sb  12-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 

>  100  knots)]  ® 

(b)  S]3  24-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 

>  100  knots)]  ^ 

(c)  S]3  36-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 

>  100  knots)]  ^ 

(d)  S|^  48-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 

>  100  knots)] 


124 


125 


126 


127 


128 


129 


130 


131 


132 


133 


134 


135 


136 


X 


Figure 

Title 

Page 

11  (a) 

Sg  12-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 
>  100  knotsij 

137 

(b) 

Sg^  24-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 
>  100  knots)]  ^ 

138 

(c) 

Sg  36-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 
>  100  knots)] 

139 

(d) 

Sg  48-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 
>  100  knots)] 

140 

12  (a) 

Sg  12-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 
>  100  knots'll  ^ 

141 

(b) 

Sg  24-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 
>  100  knots »[  ^ 

142 

(c) 

Sg  36-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 
>  100  knots)]  ^ 

143 

(d) 

Sg  48-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as 
>  100  knots)]  ^ 

144 

LIST  OF  TABLES 

Table 

Title 

Page 

I 

Grid  point  categorization 

74 

II 

Possible  predictands 

76 

III 

Predictors  available  for  screening 

78 

IV 

Reduction  of  predictands 

80 

V 

Reduction  of  predictors 

81 

VI 

Selected  predictors  for  polar  jet-stream  categories 

83 

vn 

Selected  predictors  for  subtropical  jet-stream  categories 

87 

xi 


Table  Title  Page 

VIII  Rms  errors  on  independent  data  for  polar  jet-stream  categories  90 

IX  Rms  errors  on  independent  data  for  subtropical  jet-stream 

categories  94 

X  Rms  errors  on  independent  data  for  comparison  of  category 

equations  97 

XI  Overall  rms  errors  for  LMW  prognoses  108 


PART  I 


A  TECHNIQUE  for  ANALYSIS  of 
LEVEL  of  MAXIMUM  WIND  PARAMETERS 


by 


David  B.  Spiegler 


SECTION  I 


INTRODUCTION 

The  National  Meteorological  Center  (NMC)  presently  sends  hand-analyzed  level- 
of-maximum  wind  (LMW)  analyses  for  the  United  States  and  immediate  areas  outside 
its  borders  to  users  via  the  national  facsimile  network.  These  analyses  are  for  the 
maximum  wind  speed,  the  level  of  maximum  wind,  and  an  averaged  vertical  shear  below 
and  above  the  LMW. 

The  work  described  here  represents  an  attempt  to  develop  an  objective  (computer) 
technique  for  the  analysis  of  five  LMW  parameters  for  the  Northern  Hemisphere.  These 
parameters  are: 

(a)  the  wind  direction  at  the  LMW  [W^(L)], 

(b)  the  wind  speed  at  the  LMW  [W  (L)], 

s 

(c)  the  height  of  the  LMW  [Z(L)], 

(d)  the  shear  below  the  LMW 

(e)  the  shear  above  the  LMW  (S  ). 

3. 

The  objective  technique  was  developed  in  direct  support  of  a  3rd  Weather  Wing 
requirement  for  a  program  dealing  with  an  improved  method  for  analyzing  and  fore¬ 
casting  clear-air  turbulence. 

Because  the  analyses  are  for  the  Northern  Hemisphere— an  area  that  has  widely- 
variable  data  density— the  initial-guess  fields  are  of  prime  importance.  Thus,  the  analy¬ 
sis  technique  uses  an  established  relationship^  between  the  LMW  parameters  and  the 
information  routinely  available  from  constant-pressure-surface  analysis  of  height, 
temperature,  and  wind.  The  initial  guesses  are  corrected  by  using  station  observations 
of  the  LMW  parameters  (from  a  preprocessed  data  tape)  and  the  successive  approxima¬ 
tion  analysis  technique  (SAT)  [1]. 

For  the  wind-speed  analysis  at  the  LMW  the  jet  cores  are  located,  and  “observa¬ 
tions”  are  “generated”  at  the  jet  core  by  using  models  of  horizontal  jet  profiles. 

Verifications  of  the  analyses  are  determined  by  the  areal-mean-error  method  [6]. 


The  relationship  is  established  by  means  of  regression  equations,  for  each  of  the 
LMW  parameters,  that  are  a  function  of  the  type  of  jet  and  its  strength  (e.g.,  strong- 
polar  jet,  weak -moderate  subtropical  jet,  etc.). 
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SECTION  II 


LMW  SCREENING  REGRESSION  EXPERIMENTS 

Three  data  preprocessor  programs  were  necessary  for  the  derivation  of  regres¬ 
sion  equations  that  are  used,  with  constant-pressure -surf ace  analysis,  to  generate  initial- 
guess  fields  for  the  LMW  parameters.  Descriptions  of  the  functions  of  the  three  programs 
follow: 

(a)  Jet  Stream  Preprocessor  I  (JSP-I )— combines  the  wind  data  from  the 

constant-pressure  radiosonde  observations  with  the  rawinsonde  data.  It  also 

determines  the  height  of  LMW  and  the  wind  speed  at  the  LMW,  and  computes  the 

vertical  wind  shear  above  and  below  the  LMW  for  each  of  the  observation  stations. 

If  there  are  wind  speeds  in  the  sounding  that  are  within  three  knots  of  W  (L),  the 

s 

wind  profile  is  considered  to  contain  a  multiple  maximum.  In  these  cases, 

Z(L)  is  determined  by  averaging  the  Z(L)  reports  within  one  grid  interval  of  the 
‘^problem’^  station;  the  calculated  Z(L)  at  that  station  is  the  height  that  is  closest 
to  the  average  of  the  Z(L)’s  from  the  nearby  stations. 

(b)  Jet  Stream  Preprocessor  II  (JSP-II)— produces  a  tape  with  a  total  of  58 
parameters  that  are  obtained  directly,  or  computed,  from  information  available  at 
constant  pressure  surfaces  and  which  are  considered  as  predictors  for  LMW  fields. 
These  58  parameters,  including  four  LMW  parameters  [W^(L)  is  excluded],  are: 


(i)  Z(7) 

Z(5)  Z(4)  Z(3) 

Z(2)  Z(l) 

T(7) 

T(5)  T(4)  T(3) 

T(2)  T(l) 

Wd(7) 

W^(5)  W^(4)  W^(3) 

W^(2)  W^(l) 

W  (7) 

W  (5)  W 

(4)  W  (3) 

W  (2)  W  (1) 

s 

s  s 

s 

s  s 

where:  Z  =  height;  T  = 

temperature; 

W(j  =  wind  direction;  Wg  =  wind  speed;  and  7,  5, 

4,  3,  2,  and  1  represent  the  700-,  500-, 

,  400-,  300-, 

200-,  and  100-mb  levels  [e.g.,  Z(7) 

is  the  height  of  the  700-mb  level]. 

(ii)  h(7-5) 

AT(7— 5) 

AWg(7-5) 

AT/An(7— 5) 

(9h*/9n  sin  Q!)  (5) 

h(5-4) 

AT(5-4) 

AWs(5-4) 

AT/An(5— 4) 

(9h*/9n  sin  oi)  (4) 

h(5-3) 

AT(5-3) 

AWs(5-3) 

AT/An(5-3) 

(9h*/9n  sin  oi)  (3) 

h(4-3) 

AT(4— 3) 

AWs(4-3) 

AT/An(4-3) 

(9h*/9n  sin  a)  (2) 

h(3— 2) 

AT(3-2) 

AWg(3-2) 

AT/An(3— 2) 

(9h*/9n  sin  o')  (i) 

h(2-l) 

AT(2-1) 

AWg(2-l) 

AT/An(2-1) 

2 


where:  h  =  thickness  of  a  layer  bounded  by  two  mandatory  constant  pressure 

surfaces;  AT  =  difference  in  temperature  between  two  mandatory  constant 

pressure  surfaces;  AW  =  difference  in  magnitude  of  the  wind  speed  between 

s 

constant  pressure  surfaces;  AT/ An  =  the  horizontal  temperature  gradient  in 
the  middle  of  a  layer  between  constant  pressure  surfaces;  and  9h*/9n  sin  oi  =  the 
temperature  gradients  at  a  constant  pressure  surface  where  h*  is  a  pseudo¬ 
thickness  and  sin  ot  denotes  the  direction  of  the  gradient  {ot  -  the  angle  between 
the  thermal  and  actual  wind). 

(ill)  Z(L)  W^(L) 

where  r  =  the  radius  of  curvature  at  200  mb.  The  purpose  of  this  program  is 
to  have  the  data  in  a  form  that  can  be  used  with  a  screening  regression  program 
to  generate  equations  for  specifying  the  four  LMW  parameters  [W^(L)  is  speci¬ 
fied  separately  by  a  procedure  described  in  Section  III]. 

(c)  Jet  Stream  Categories  (JSC)— stratifies  the  station  data  on  the  basis 
of  wind  speeds  at  200  and  300  mb,  the  latitude  of  the  station,  and  a  combination 
of  both  wind  speeds  and  latitude.  The  purpose  of  this  program  is  to  have  the 
screening  regression  derive  different  equations  for  jets  of  various  strength  and 
type.  The  exact  specification  of  the  categories  (see  Table  I)  was  determined  by 
examining  several  months  of  LMW  data. 

Three  months  of  a  winter  upper-air  data  sample  from  the  National  Meteorological 
Center  were  run  through  the  three  preprocessor  programs.  Screening  regression  experi¬ 
ments  for  obtaining  equations  for  the  initial  guesses  of  LMW  parameters  were  designed. 
The  main  purpose  of  the  experiments  was  to  obtain  the  most  practical  set  of  equations 
from  an  operational  standpoint,  while  keeping  the  errors  within  acceptable  limits. 
Selective  “addition”  and  “subtraction”  of  possible  predictors  [listed  in  (b)  above]  for 
several  different  screening  runs  gave  the  information  necessary  to  decide  which 
equations  were  the  most  useful. 

A  program  that  derives  multiple  linear-regression  equations  by  using  a  screening 
procedure  has  been  written  and  was  used  in  the  study.  The  screening  procedure  selects, 
from  a  set  of  possible  predictors,  a  subset  that  contributes  significantly  and  independ¬ 
ently  to  reducing  the  variance  of  the  predictand.  The  predictor  having  the  highest  linear 
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correlation  with  the  specified  predictand  is  selected  from  an  array  of  possible  pre¬ 
dictors.  Next,  the  partial  correlation  coefficients  between  each  of  the  remaining 
predictors  and  the  predictand  (holding  the  first  predictor  constant)  are  examined,  and 
the  predictor  associated  with  the  highest  coefficient  is  then  selected  as  the  second 
predictor.  Additional  predictors  are  selected  in  a  similar  manner.  This  procedure  is 
repeated  until  a  chosen  predictor  fails  to  explain  a  significant  additional  percentage  of 
the  remaining  variance  of  the  predictand  or  until  a  specified  number  of  predictors  is 
selected.  The  criterion  of  significance  employed  is  a  modified  F-test  described  by 
Miller  [3]. 

To  justify  the  definitions  used  for  the  various  LMW  categories,  contingency  tables 
were  computed  that  relate  the  category  definitions  to  definitions  of  jet  types  gained 
by  synoptic  experience  and  given  by  the  height  of  the  LMW. 

For  the  purposes  of  the  contingency  tables,  the  definitions  of  polar  and  sub¬ 
tropical  jets,  and  jets  that  are  not  clearly  polar  or  subtropical  (subtropical-polar)  are 
as  follows  for  winter  data: 

(a)  Polar  [Z(L)  <  34,000  ft] 

(b)  Subtropical  [Z(L)  >  37,000  ft] 

(c)  Subtropical— polar  [34,000  Z(L)  :£  37,000  ft] 

The  contingency  tables  [Table  II,  (a)  and  (b)]  are  based  on  station  data  from  the 
developmental  data  sample  (December  1963,  and  January— February  1964).  The  defini¬ 
tions  based  on  300-  and  200-mb  wind  speeds  and  latitude  appear  to  be  satisfactory  with 
only  a  small  percentage  of  the  cases  in  the  lower  left  or  upper  right  corners  of  the 
tables. 
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SECTION  m 


METHOD  OF  ANALYSIS  AND  VERIFICATIONS  OF  LMW  PARAMETERS 

Initial-guess  fields  for  all  LMW  parameters  except  W^(L),  are  generated  by  the 

derived  regression  equations.  The  analysis  program  determines  which  jet-stream 

category  a  particular  grid  point  is  in,  and  uses  the  appropriate  regression  equations  for 

each  of  the  LMW  parameters.  The  initial  guesses  are  not  generated  simultaneously, 

because  one  (or  more)  already-analyzed  LMW  fields  are  often  predictors  in  the  equation(s) 

for  the  other  field(s).  The  initial  guesses  are  adjusted  by  using  the  preprocessed  LMW 

station  data  and  the  successive  approximation  analysis  technique  [1].  For  the  W  (L) 

s 

field,  the  jet  core  is  located,  and  observations  are  generated  at  the  core  by  employing 
models  of  horizontal  jet-stream  profiles.  The  details  for  the  determination  of  each  of 
the  LMW  fields,  their  order  in  the  analysis  procedure,  and  the  verification  procedure 
follow. 

1.  Wind  Direction  at  the  LMW  [W,(L)] 

d 

The  procedure  for  determining  wind  direction  at  the  LMW  is: 

(a)  for  grid  points  that  are  at  latitude  <  30°N,  the  wind  direction  (W^)  at 

200  mb  [W,(2)]  is  used  as  the  W,(L). 
d  d 

(b)  for  grid  points  that  lie  between  30°  and  45°N:  if  the  radius  of  curvature 
at  200  mb  (r  )  is  positive,  W  (3)  is  used  as  W  (L);  if  r  is  negative  or  zero, 

^  Q  G.  Z 

W^(2)  is  specified  as  the  W^(L). 

(c)  for  grid  points  that  are  at  latitudes  >  45°N,  W^(3)  is  also  the  W^(L). 

2.  Wind  Speed  at  the  LMW  [W  (L)] 

s 

The  initial-guess  field  is  generated  by  the  appropriate  equations  (see  the  Appendix) 

that  are  dependent  upon  the  category  in  which  a  particular  grid  point  is  located.  The 

initial-guess  field  of  W  (L)  is  “corrected”  by  using  station  observations  of  W  (L)  and 

s  s 

SAT.  The  number  of  passes  through  the  data  is  an  input  item,  but  the  maximum  number 
of  passes  is  five,  some  of  which  may  be  made  with  additional  generated  observations. 

The  jet  cores  are  located  after  a  specified  number  of  passes.  The  first  approxi¬ 
mation  to  defining  the  jet  core  is  obtained  by  searching  among  grid  points  for  the  maxi- 
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mum  wind  speed  within  a  specified  area.  For  the  analyses  using  the  winter-data  sample 
the  areas  are  enclosed  by  those  grid  points  having  a  wind-speed  value  greater  than  or 
equal  to  75  knots. 

A  systematic  search  is  performed  within  the  area  by  tagging  grid  points  that  are 
maximums  (the  search  is  systematic  in  the  sense  that  it  proceeds  independently  along 
the  vertical,  horizontal,  and  diagonal  lines  of  grid  points).  The  grid  points  with  the  most 
“tags”  are  the  first  approximations  to  the  jet  cores.  More  specific  details  and  a  discus¬ 
sion  of,  and  solutions  for,  the  various  problems  are  contained  in  the  specifications  for 
the  LiMW  computer  program  [5]. 

The  set  of  consecutive  stored  grid  points,  which  are  the  first  approximation  to  a 
jet  core,  is  then  used  with  a  parabola-fitting  routine  to  determine  the  “exact”  location 
of  a  jet  core  between  grid  points. 

After  the  jet  cores  have  been  located,  observations  are  generated  at  the  midpoints 
of  the  jet-core  segments  by  using  Endlich  and  McLean’s  jet-stream  model  [2]  which 
relates  the  wind  speed  at  the  core  to  the  wind  speed  at  some  specified  perpendicular 
distance  from  the  core  (see  Fig.  1).  By  using  the  model,  two  observations  are  generated 
at  the  same  points  on  the  core,  i.e.,  one  from  the  wind  on  the  high-pressure  side  and  the 
other  from  the  wind  on  the  low-pressure  side.  The  average  of  these  two  “generated” 
observations  is  used  as  the  final  observation  for  the  midpoint  of  a  core  segment. 

There  is  an  option  to  make  an  additional  pass  for  analysis  of  W  (L)  with  the 

s 

generated  data  and  the  actual  station  data  using  a  relatively-small  influence 

radius.  The  generated  observations  may  be  assigned  a  weight  relative  to  the  actual  data 

for  this  pass. 

3.  Verification  of  W  (L) 

s  ' 

The  analysis  of  W  (L)  is  verified  by  the  areal-mean/ root-mean-square  error 
s 

method  [6].  This  method  approximates  an  area  integral  of  the  error  over  the  analysis 
area  by  withholding  some  percentage  of  the  observations  and  verifying  them  against  the 
withheld-  and  analysis-station  observations,  which  are  equally  weighted.  The  withhe Id- 
station  observations  approximate  the  maximum  error  points  in  the  analysis  field,  but 
the  number  of  withheld  stations  required  by  the  areal-mean-error  method  usually  does 
not  equal  the  number  of  maximum-error  points.  Therefore,  the  withheld  stations  were 
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Fig.  1.  Average  horizontal  profile  of  wind  speeds  in  the  jet  stream  (taken  from  [2]). 


selected  to  represent  equally  the  entire  spectrum  of  maximum  errors  over  the  analysis 
area.  The  assumption  is  that  the  maximum  error  is  inversely  proportional  to  the  density 
of  the  data.  The  density  of  each  station  is  computed  prior  to  the  analysis,  and  withheld 
stations  are  selected  on  the  basis  of  three  LMW-wind-speed  and  three  station-density 
categories  (see  Table  III). 


TABLE  III 

SELECTION  OF  WITHHELD  STATIONS  ON  THE  BASIS  OF 
STATION  DENSITY  AND  LMW  WIND  SPEED 


Density  (p^) 

Category 

A 

No.  of  stations 

W  (L)  >  100 

B 

No.  of  stations 

51  ^  W  (L)  ^  100 

C 

No.  of  stations 

W  (L)  50 

s 

Low;  0.11  to  0.32 

3 

4 

3 

Med;  0.33  to  0.43 

4 

3 

3 

High;  s  0.44 

3 

3 

4 

4.  Analysis  and  Verification  of  Z(L) 

The  initial-guess  field  for  Z(L)  is  generated  with  the  appropriate  regression  equa¬ 
tions  after  determining  which  category  each  grid  point  is  in.  The  successive  approxima¬ 
tion  analysis  technique  is  then  used  with  the  station  data  for  Z(L)  to  adjust  the  initial- 
guess  field. 

Where  the  W  (L)  is  less  than  75  knots,  Z(L)  at  adjacent  grid  points  and/or 
s 

stations  may  vary  widely;  this  is  a  common  characteristic  of  adjacent  wind  profiles  that 

exhibit  relatively  weak  maximums.  For  this  reason,  the  analysis  is  meaningful  only 

over  those  areas  on  the  Z(L)  chart  where  the  W  (L)  is  greater  than  75  knots. 

s 

Verification  of  Z(L)  is  done  with  the  areal-mean/root-mean-square  error  method 

using  the  identical  set  of  withheld  stations  used  in  the  verification  of  W  (L). 

s 

5.  Analysis  and  Verification  of  and 

Both  of  these  parameters  have  initial  guesses  generated  by  the  appropriate  regression 

equations.  As  may  be  seen  from  the  equations  in  the  Appendix,  W  (L)  and  Z(L)  are  often 

s 

“predictors”  for  S|^  and  Sg. 

Once  again,  SAT  is  used  with  the  preprocessed  station  data  to  adjust  the  initial- 
guess  fields. 
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SECTION  IV 


RESULTS 

The  evaluation  of  the  quality  of  the  analysis  technique  was  based  on  the  results 
from  two  general  classes  of  experiments— the  experiments  for  deriving  the  initial- 
guess  equations  for  the  LMW  analysis  fields,  and  the  experiments  with  the  analyses 
that  were  generated  by  the  analysis  procedure. 

6.  Results  of  LMW  Regression  Experiments 

Non-derived  predictors  (e.g.,  temperatures  and  wind  speeds)  and  derived  predictors 

2 

(e.g.,  thermal  winds  and  shears)  are  among  the  possible  predictors  for  the  LMW  param¬ 
eters  (see  Section  I).  In  the  design  of  the  experiments,  a  minimum  of  three  to  a  maximum 
of  seven  screening  runs  were  specified  for  each  predictand  in  each  category.  The  dif¬ 
ferences  between  the  individual  screening  runs  for  a  particular  predictand  were  in  the 
variables  specified  as  possible  predictors.  The  objective  in  having  several  screening 
runs  was  to  answer  the  following  questions: 

(a)  Will  the  non-derived  (simple)  predictors  alone  satisfactorily  specify 
the  LMW  parameters? 

(b)  Will  the  derived  (sophisticated)  predictors  alone  give  a  better  result 
than  the  non-derived  predictors  alone? 

(c)  Will  a  combination  of  non-derived  and  derived  predictors  yield  the 
best  results? 

(d)  Does  the  height  at  the  LMW  [Z(L)]  add  more  to  the  specification  of  the 

wind  speed  at  the  LMW  [W  (L)]  than  the  wind  speed  adds  to  the  height? 

s 

(e)  Are  the  previously  analyzed  LMW  fields  important  in  specifying  the 
shears  above  and  below  the  LMW? 

The  answers  to  these  questions  are  contained  in  Table  IV,  which  lists  the  selected 
predictors  for  the  LMW  parameters  for  each  of  the  LMW  categories,  and  the  associated 
dependent  and  independent  data  results.  The  LMW  screening  regression  equations  for 
each  type  of  jet  category  are  in  the  Appendix. 


W^  and  Z  at  constant-pressure  surfaces,  and  all  parameters  containing  400-mb 
information,  were  not  included  as  possible  predictors  in  any  of  the  screening  runs. 
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TABLE  IV 

INITIAL-GUESS  EQUATIONS  FOR  LMW  PARAMETERS- 
DEPENDENT  AND  INDEPENDENT  DATA  RESULTS  FOR  SEVEN  JET  STREAM  CATEGORIES 


Category 

Predlctand 

Selected  predictors 

%  red. 

Number 
of  cases 

rms  error 

Std.  dev. 

Means 

Units 

(dep) 

(dep) 

(Ind) 

(dep) 

(Ind) 

(dep) 

(Ind) 

(dep) 

(Ind) 

Z(L) 

Wg(l).  W^(7).  T(7) 

37 

3482 

467 

9611 

9472 

12076 

12630 

34403 

33850 

ft 

W,(L) 

W^(2),  W^(l),  W^(5) 

60 

3482 

467 

7.5 

9.7 

11.9 

12.2 

42.7 

43.9 

knot 

1 

s 

a 

Sb,  W^(l),  W^(L) 

39 

3482 

467 

0.93 

0.95 

1.19 

1.20 

-1.48 

-1.87 

knot  10“^  ft 

W^(L).  W^(2),  W^(5),  Z(L) 

41 

3482 

467 

0.79 

1.48 

1.04 

1.46 

+1.72 

+1.79 

knot  10"^  ft 

Z(L) 

Wg(l),  W^(3),  T(5) 

14 

2010 

317 

4740 

3806 

5127 

4199 

40019 

40125 

ft 

W^(L) 

W^(2) 

55 

2010 

317 

12.5 

8.8 

18.6 

16.0 

83.5 

82.6 

knot 

2 

s 

a 

W  (L),  W  (1),  W  (2) 
s  s  s 

32 

2010 

317 

2.49 

1.29 

3.01 

1.61 

-3.14 

-2.73 

knot  10  ^  ft 

Z(L).  W^(L),  W^(3) 

63 

2010 

317 

2.14 

1.18 

3.53 

1.58 

+3 .32 

+3.01 

knot  10  ^  ft 

Z(L) 

Wg(3),  W^(l),  W^(2) 

06 

1141 

135 

3914 

3291 

4046 

3596 

39214 

38750 

ft 

W^(L) 

W^(2) 

67 

1141 

135 

12.7 

10.0 

22.0 

24.2 

132.7 

132.07 

knot 

3 

s 

a 

W^(L).  W^(2),  Z(L).  Wg{l) 

39 

1141 

135 

2.98 

1.83 

3.83 

2.75 

-4.96 

-4.95 

knot  10*^  ft 

S 

a 

W^{L),  W^(2),  Z(L),  W^(l). 

43 

1141 

135 

2.88 

1.72 

3.83 

2.75 

-4.96 

-4 .95 

knot  10  ^  ft 

Z(L),  W^(3),  W^(L) 

81 

1141 

135 

1.24 

1.13 

2.84 

2.38 

+4.57 

+4.28 

knot  10  ^  ft 

Z(L) 

W^(l),  T(7).  W^(5),  W^(3),  W^(2) 

23 

1363 

211 

5767 

5512 

6587 

6896 

35659 

37584 

ft 

W^(L) 

W^(2) 

54 

1363 

211 

12.7 

8.9 

18.6 

15.9 

80.5 

81.8 

knot 

4 

S 

a 

W^(L),  W^(3),  W^(l) 

36 

1363 

211 

2.34 

1.33 

2.94 

1.34 

-2.63 

-2.32 

knot  10  ^  ft 

\ 

W^(L),  Z(L),  W^{3),  W^(5) 

49 

1363 

211 

2.50 

1.14 

3.48 

1.39 

+2.55 

+2.24 

knot  lo’^  ft 

Z(L) 

T(7),  W^(l).  W^(3),  W^(2) 

19 

477 

70 

4736 

4668 

5251 

5061 

36339 

35494 

ft 

W^(L) 

W^(3),  W^(2) 

71 

477 

70 

11.2 

10.2 

21.0 

22.6 

132.2 

126.84 

knot 

5 

S 

a 

Z(L),  S^,  W^(l) 

50 

477 

70 

1.87 

2.04 

2.66 

2.25 

-3.74 

-3.22 

knot  10  ^  ft 

S 

a 

Z(L),  S^,  W^(l),  W^(L),  W^(2) 

5B 

477 

70 

1.73 

1.80 

2.66 

2.25 

-3.74 

-3.22 

knot  lO"^  ft 

Sb 

Z(L).  W^(L),  W^(3) 

69 

477 

70 

1.59 

1.42 

2.87 

2.66 

+3.57 

+3.57 

knot  10  ^  ft 

Z(L) 

Wg(5).  W^{2).  W^(3) 

27 

2543 

319 

4991 

6440 

5830 

7396 

30996 

33459 

ft 

W^(L) 

W^{3),  W^(2) 

68 

2543 

319 

8.7 

11.3 

15.4 

16.6 

77.4 

80.8 

knot 

6 

S 

a 

W^{L),  W^(2).  Z(L).  W^(l) 

56 

2543 

319 

1.09 

1.46 

1.65 

2.06 

-3.06 

-3.07 

knot  lO"^  ft 

Wg(L),  W^<5),  Z(L).  W^(3) 

67 

2543 

319 

0.93 

0.91 

1.61 

1.63 

+2.97 

+2.97 

knot  10  ^  ft 

Z(L) 

T<5),  W^<2).  W^<3) 

26 

477 

56 

3407 

4429 

3958 

4909 

30684 

31701 

ft 

W3(L) 

W^(3),  W^(2) 

64 

477 

56 

11.6 

12.3 

19.4 

22.5 

126.5 

130.4 

knot 

7 

s 

a 

W^(L),  Wg(2),  Z(L),  Wg(3) 

67 

477 

56 

1.64 

2.00 

2.87 

2.43 

-5.55 

-4.90 

_3 

knot  10  ft 

Z(L).  W^{L),  W^(5) 

62 

477 

56 

1.68 

1.90 

2.73 

2.30 

+5.16 

+5.14 

knot  10  ^  ft 

- 

Z(L),  W^(L),  W^(5).  W^(3) 

67 

477 

56 

1.57 

1.72 

2.73 

2.30 

+5.16 

+5.14 

knot  10  ^  ft 
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The  outcome  of  the  experiments,  with  respect  to  the  selected  predictors,  is 
interesting.  Percent  reduction  of  variance  {%  red.)  and  rms  errors  were  the  criteria 
used  for  choosing  the  equations.  It  is  evident  from  Table  IV  that  the  derived  predictors 
failed  to  add  significantly  to  the  non-derived  predictors.  It  is  also  evident  that  the 
derived  predictors  alone  did  not  yield  results  as  good  as  those  from  the  non-derived 
predictors.  Table  V  indicates  the  superiority  of  the  non-derived  predictors  for  each 
predictand  and  each  category. 

For  jet-stream  categories  2  through  7,  the  screening  runs  using  non-derived 
predictors,  from  which  the  equations  were  selected,  are  compared  with  screening  runs 
that  use  derived  predictors  (it  is  implicit  that  if  LMW  predictors  were  used  in  the 
experiment  for  the  simple  predictors,  the  same  LMW  predictors  were  used  in  the  experi¬ 
ment  with  the  sophisticated  predictors). 

A  smaller  number  of  experiments  were  run  for  Category  1  than  for  Categories  2 
through  7.  Table  V  contains  all  the  experiments  for  the  seven  categories.  There  is  no 
direct  comparison  of  predictors  in  Category  1,  but  results  indicated  that  non-derived 
predictors  were  more  highly  correlated  with  the  LMW  parameters  than  were  the  derived 
predictors. 

The  following  conclusions  may  be  drawn  from  the  table: 

(a)  A  higher  number  of  derived  predictors  is  often  necessary  to  attain  the 
same  percent  reduction  as  the  simple  predictors. 

(b)  For  two  thirds  of  the  predictands  (Categories  2—7),  an  the  derived 
predictors  (plus  LMW  parameters  where  indicated)  do  not  result  in  as  high  a 
percent  reduction  as  four  or  less  simple  predictors  (plus  LMW  parameters  where 
indicated). 

(c)  When  the  number  of  derived  and  non-derived  predictors  is  the  same,  the 
percent  reductions  are  within  2.5%  of  each  other.  In  Categories  2—7  this  occurs 
only  with  Z(L).  Although  the  percent  reductions  are  approximately  the  same  for 
these  six  cases,  it  is  easier  to  use  non-derived  fields  in  an  operational  system. 

The  results  for  the  LMW  equations  on  independent  data  (Table  IV)  indicate  that 

the  equations  are  stable,  with  independent  errors  near  or  lower  than  the  dependent 
errors. 
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TABLE  V 

PERCENT  REDUCTION  OF  VARIANCE  ASSOCIATED  WITH 
DERIVED  AND  NON-DERIVED  PREDICTORS* 

(a)  Category  1 


Predictand 

Possible 

Selected 

% 

Total 

predictors 

predictors 

red. 

%  red. 

Ws(l) 

19.6 

NDt 

W^(7) 

12.1 

36.7 

T(7) 

5.0 

Z(L) 

Ws(l) 

19.6 

D  +  ND 

W^(7) 

12.1 

37.5 

AW  (3-2) 
s 

5.8 

Ws(2) 

40.0 

NDt 

Ws(l) 

12.7 

60.4 

W  (L) 
s  ' 

Ws(5) 

7.7 

W^(2) 

40.0 

D  +  ND 

AW  (2—1) 

12.7 

60.4 

Ws(5) 

7.7 

W^(l) 

3.1 

ND 

W^(3) 

4.1 

8.5 

T(3) 

1.3 

AWg(2— 1) 

3.2 

S 

a 

D  +  ND 

AW  (3-2) 
s' 

4.0 

8.3 

W^(7) 

1.1 

12.8 

D  +  ND  +  LMWt 

W^(l) 

5.0 

38.7 

Ws(L) 

20.9 

Ws(l) 

1.4 

ND 

T(5) 

1.1 

2.5 

T(l) 

0.0 

Ws(l) 

1.4 

D  +  ND 

AW  (5-3) 
W^(3) 

0.2 

0.0 

1.6 

W  (L) 

20.2 

D  +  ND  +  LMWt 

W^(2) 

8.1 

41.1 

Ws(5) 

5.9 

Z(L) 

6.9 

*D  =  derived;  ND  =  non-derived. 
tSelected  for  use. 
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TABLE  V  (cont’d) 


(b)  Category  2 


Predictand 

Possible 

Selected 

% 

Total 

predictors 

predictors 

red. 

%  red. 

W^(l) 

6.1 

NDt 

W^(3) 

5.4 

14.5 

Z(L) 

T(5) 

3.0 

W^(2-l) 

7.2 

AW^(3-2) 

4.2 

D 

AT(3-2) 

1.1 

14.8 

AW  (2-1) 

1.7 

W^(5-3) 

0.6 

W^(L) 

NDt 

Ws(2) 

55.1 

55.1 

D 

ALL 

28.9 

28.9 

Ws(L) 

21.3 

ND  +  LMWt 

Ws(l) 

7.5 

31.8 

W^(2) 

3.0 

s 

a 

W^(L) 

21.3 

D  +  LMW 

^b 

4.5 

33.3 

Z(L) 

4.6 

AW  (2—1) 

S 

2.9 

Z(L) 

28.6 

ND  +  LMWt 

Ws(L) 

15.2 

63.5 

W^(3) 

19.7 

D  +  LMW 

ALL 

58.5 

58.5 
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TABLE  V  (cont’d) 


(c)  Category  3 


Predictand 

Possible 

Selected 

% 

Total 

predictors 

predictors 

red. 

%  red. 

Ws(3) 

2.1 

NDt 

W^(l) 

3.1 

6.4 

Z(L) 

Ws(2) 

1.2 

AW^(3-2) 

4.4 

D 

W^(2-l) 

2.9 

8.3 

h(7-5) 

1.0 

W  (L) 

NDt 

Ws(2) 

66.5 

66.5 

s 

D 

ALL 

49.5 

49.5 

Ws(L) 

14.3 

W  (2) 

11.4 

S 

a 

ND  +  LMWt 

s 

Z(L) 

6.5 

39.2 

W^(l) 

7.0 

W^(L) 

6.0 

Z(L) 

12.7 

D  +  LMW 

3.8 

40.3 

AW  (3—2) 
s 

3.5 

AWg(2-l) 

14.3 

Z(L) 

25.8 

ND  +  LMWt 

Ws(3) 

W^(L) 

11.6 

43.6 

81.0 

D  +  LMW 

all 

71.3 

71.3 
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TABLE  V  (cont’d) 


(d)  Category  4 


Predictand 

Possible 

Selected 

% 

Total 

predictors 

predictors 

red. 

%  red. 

W^(l) 

4.8 

T(7) 

5.2 

NDt 

W^(5) 

4.6 

23.3 

W^(3) 

1.4 

W  (2) 

7.3 

Z(L) 

s 

AW^3-2) 

13.8 

W^(2-l) 

3.4 

D 

W^(5-3) 

3.4 

25.1 

Ah(7-5) 

3.9 

AT(2— 1) 

0.6 

W  (L) 

NDt 

W  (2) 

53.6 

53.6 

s 

D 

ALL 

28.9 

28.9 

W^(L) 

22.5 

S 

ND  +  LMWt 

W^(3) 

10.9 

36.5 

a 

W^d) 

3.1 

D  +  LMW 

ALL 

34.5 

34.5 

W^(L) 

17.9 

ND  +  LMWt 

Z(L) 

14.7 

48.6 

W^(3) 

11.9 

W^(2) 

4.1 

D  +  LMW 

ALL 

40.0 

40.0 
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TABLE  V  (cont’d) 


(e)  Category  5 


Predictand 

Possible 

Selected 

% 

Total 

predictors 

predictors 

red. 

%  red. 

T(7) 

4.0 

Ws(l) 

2.7 

NDt 

Ws(3) 

3.8 

18.6 

Z(L) 

Ws(2) 

8.1 

AW  (3-2) 
s 

13,5 

D 

W^(2-l) 

3.4 

19.4 

Ah(5-3) 

1.7 

"2 

0.8 

Ws(3) 

68.2 

W^(L) 

NDt 

W^(2) 

3.3 

71.5 

D 

ALL 

50.6 

50.6 

Z(L) 

12.6 

s. 

29.4 

s 

ND  +  LMWt 

D 

W^(l) 

8.2 

57.6 

a 

W^(L) 

2.0 

W^(2) 

5.4 

D  +  LMW 

ALL 

57.2 

57.2 

Z(L) 

35.9 

S. 

ND  +  LMWt 

W^(L) 

15.8 

69.2 

D 

W^(3) 

17.5 

D  +  LMW 

ALL 

63.7 

63.7 
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TABLE  V  (cont’d) 


(f)  Category  6 


Predictand 

Possible 

Selected 

% 

Total 

predictors 

predictors 

red. 

%  red. 

W^(5) 

7.5 

NDt 

W^(2) 

13.6 

26.7 

Z(L) 

Ws(3) 

5.6 

W^(3-2) 

22.7 

D 

W^(5-3) 

4.0 

29.1 

AT(5-3) 

2.4 

W^(3) 

59.9 

W  (L) 
s 

NDt 

8.2 

68.1 

D 

ALL 

35.1 

35.1 

W^(L) 

22.1 

ND  +  LMWt 

Ws(2) 

21.2 

55.9 

s 

Z(L) 

5.4 

a 

W^(l) 

7.2 

D  +  LMW 

ALL 

46.1 

46.1 

W^(L) 

21.2 

ND  +  LMWt 

Ws(5) 

20.1 

66.8 

Z(L) 

20.6 

W  (3) 

4.9 

D  +  LMW 

ALL 

49.0 

49.0 
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TABLE  V  (cont’d) 


(g)  Category  7 


Predictand 

Possible 

predictors 

Selected 

predictors 

% 

red. 

Total 
%  red. 

T(5) 

14.6 

NDt 

Ws(2) 

5.3 

25.9 

Z(L) 

W^(3) 

6.0 

W^(3-2) 

17.9 

D 

Ah  (7 -5) 

6.3 

25.9 

AT(3-2) 

1.7 

W^(L) 

NDt 

Ws(3) 

W^(2) 

57.7 

6.9 

64.6 

D 

ALL 

48.2 

48.2 

W^(L) 

15.1 

ND  +  LMWt 

W^(2) 

25.5 

67.4 

s 

a 

Z(L) 

20.6 

W^(3) 

6.2 

D  +  LMW 

ALL 

64.8 

64.8 

Z(L) 

16.5 

S. 

ND  +  LMWt 

W^(L) 

16.0 

62.3 

D 

W^(5) 

29.8 

D  +  LMW 

ALL 

52.6 

52.6 
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The  above  discussion  of  the  tables  answers  the  first  three  questions  posed 

earlier  in  this  section.  Answers  to  questions  (d)  and  (e)  are  found  in  Table  IV  under 

the  column  headed  “Selected  predictors”. 

The  height  of  the  LMW  is  not  selected  as  a  significant  predictor  for  W  (L),  and 

s 

vice  versa.  Both  Z(L)  and  W  (L)  are  often  selected  as  predictors  for  the  shears,  and 

s 

the  shear  below  is  selected  as  a  predictor  for  the  shear  above  in  Categories  1,  3,  and  5. 

7.  Verification  and  Evaluation  of  LMW  Analyses 

The  analysis  procedure  described  in  Section  III  was  used  with  a  winter  data 
sample  (for  a  time  period  not  within  the  three  month  sample  used  to  derive  the  equa¬ 
tions).  LMW  analyses  for  twelve  observation  times  were  generated,  verified,  and 
evaluated. 

The  verification  of  the  analyses  was  accomplished  by  using  the  areal-mean-error 
method,  and  the  evaluation  consisted  of;  (1)  a  diagnosis  of  the  meaning  of  the  error 
statistics,  (2)  a  comparison  of  the  computer  analyses  with  hand-analyzed  charts,  and 
(3)  a  consistency  check  to  determine  if  the  analyses  of  the  LM  vV  parameters  for  each 
observation  time  were  consistent  with  one  another. 

8.  Areal-mean-error -method  Statistics 

The  overall  rms  errois for  four  LMW  parameters  for  twelve  observation  times 

during  January  1963  are  presented  in  Table  VI.  The  column  headed  “Analysis  rms 

error”  shows  the  root-mean-square  error  at  stations  used  in  the  analysis,  and  the 

column  headed  “Withheld  rms  error”  shows  the  error  at  the  withheld  stations. 

The  total  rms  error  of  15.4  knots  (Category  A)  for  W  (L)  at  the  final  pass  of  the 

s 

analysis  is  quite  low  when  one  considers  that  this  category  comprises  stations 
that  are  reporting  wind  speeds  in  excess  of  100  knots.  The  error  of  the  wind  speed  for 
this  category  is  something  less  than  15%,  and  is  probably  closer  to  10%  because  there 
are  frequently  observed  maximum  winds  of  150  to  200  knots,  especially  over  Eastern 
Asia  and  the  Western  Pacific  Ocean. 

It  is  readily  evident  in  the  overall  verification  statistics  for  all  the  LMW  param¬ 
eters  that  the  reduction  in  the  total  rms  error  from  the  initial  guess  to  the  final  pass 
is  due  almost  entirely  to  the  reduction  in  the  rms  error  at  analysis  stations  (see 
Analysis-rms-error  column  of  Table  VI).  In  fact,  the  withheld-station  rms  error  is. 
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TABLE  VI 

OVERALL  rms  ERRORS  FOR  LMW  PARAMETERS  FOR 
12  OBSERVATION  TIMES 


Error  (initial  guess/final  pass) 

Parameter 

Verification 

Total 

Analysis 

Withheld 

category 

rms 

rms 

rms 

error 

error 

error 

A 

21.9 

22.7 

21.1 

15.4 

7.0 

20.6 

W  (L) 
s 

B 

12.5 

12.2 

12.8 

(knot) 

9.7 

3.5 

13.3 

10.1 

7.8 

11.9 

8.7 

2.6 

12.1 

A 

4818 

4798 

4837 

3533 

1566 

4744 

Z(L) 

TJ 

7378 

6769 

7940 

(ft) 

Xj 

5738 

2301 

7782 

9933 

10262 

9593 

7257 

3671 

9585 

A 

2.11 

1.72 

2.45 

/V 

1.75 

0.46 

2.43 

s. 

B 

1.61 

1.48 

1.73 

_3 

(knot  10  ft) 

1.22 

0.37 

1.68 

(~i 

1.05 

1.15 

0.94 

0.70 

0.31 

0.94 

A 

2.26 

1.86 

2.60 

/V 

1.87 

0.38 

2.62 

S 

B 

1.47 

1.38 

1.56 

^  -3 

(knot  10  ft) 

1.17 

0.31 

1.62 

p 

0.94 

0.77 

1.08 

0.77 

0.25 

1.05 
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in  a  few  cases,  slightly  higher  for  the  final  pass  than  it  was  for  the  initial  guess.  This 
result  is  not  surprising  because: 

(a)  The  analysis-station  rms  error  essentially  indicates  how  well  the 
analysis  technique  fits  the  observations  (and  these  points  may  be  considered 
as  minimum  error  points  over  the  analysis  area). 

(b)  The  withheld  stations,  which  approximate  maximum  error  points, 
do  not  contribute  to  the  analysis,  by  definition.  A  station  withheld  over  a 
sparse  data  area  may  be  the  only  station  in  the  area  (e.g.,  a  stationary  weather 
ship  in  the  ocean)  and  the  analysis  over  that  area  will  remain  unchanged  from 
the  initial  guess  to  the  final  pass. 

(c)  Withheld  stations  over  medium  or  high  data-density  areas  are  not 
allowed  to  influence  the  analysis.  The  surrounding  analysis  stations,  however, 
may  change  the  analysis  to  the  extent  that  the  point  at  which  there  is  a  with¬ 
held  station  observation  may  actually  be  in  less  agreement  with  the  final  analysis 
than  it  was  with  the  initial-guess  analysis. 

An  important  point  about  the  verification  statistics  is  that  prior  to  use  of  SAT 
for  analysis  with  station  data,  the  overall  rms  errors  for  the  initial  guess  of  all  LMW 
parameters  was  approximately  equal  for  both  the  analysis  stations  and  the  withheld 
stations  (approximate  minimum-  and  maximum-error  points,  respectively).  Two  inter¬ 
pretations  of  the  data  in  Table  VI  that  may  be  made  here  are:  (1)  a  good  initial  guess 
is  essential  to  any  hemispheric  analysis  procedure,  and  (2)  the  rms-error  values  for 
the  initial  guesses  of  the  LMW  parameters  indicate  high-quality  initial  guesses. 

Figure  2  clearly  illustrates,  for  a  series  of  individual  observation  times,  the 
lack  of  bias  between  the  initial-guess  Analysis-rms  and  Wthheld-rms  errors  for 

Category  A  of  W  (L).  The  Withheld-rms  error  is  significantly  lower  than  the  Analysis- 
s 

rms  error  six  times  and  the  converse  is  true  four  times;  the  remaining  two  instances 
have  insignificant  differences. 

Figure  3  shows  the  LMW  wind-speed  analyses  for  8  and  9  January,  1963.  The 
highest  initial-guess  Withheld-rms  error  occurred  on  8  January  and  the  lowest  occurred 
on  9  January.  Withheld  station  observations  are  plotted  on  the  analyses.  Figures  4  and  5 
are  the  300-  and  200-mb  height  analyses  for  8  and  9  January,  1963,  respectively.  They 
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Fig.  2.  Initial  guess  rms  errors  for  LMW  wind  speed  at  analysis  and  withheld  stations  for  12  observation  times 
(area  where  Wx  >100  knots). 


Fig,  3(a).  Objective  LMW  wind  speed  analysis,  1200Z,  January  8,  1963  (knots). 
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Fig.  3(b).  Objective  LMW  wind  speed  analysis,  1200Z,  January  9,  1963  (knots). 
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Fig,  4(a).  300-mb  analysis,  1200Z,  January  8,  1963  (10^  ft). 
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Fig.  4(b).  200-mb  analysis,  1200Z,  January  8,  1963  (10 


Fig.  5(a).  300-mb  analysis,  1200Z,  January  9,  1963  (10^  ft). 
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are  shown  because  the  predictors  for  W  (L)  for  Categories  2—7  are  the  300-  and/or 

s 

200-mb  wind  speeds  (which  are  derived  from  the  height  analyses). 

Because  vertical  wind  profiles  sometimes  indicate  high  wind  speeds  over  a 

layer  that  may  be  as  much  as  10,000  feet  or  more  thick,  with  more  than  one  identical 

maximum  speed  within  that  layer,  Z(L)  is  sometimes  difficult  to  specify.  With  this  in 

mind,  the  total  rms  error  for  verification  Category  A  [W  (L)  >100  knots]  for  3533  ft 

s 

is  considered  to  be  a  reasonably  low  value. 

Where  W  (L)  is  less  than  100  knots  and,  particularly,  less  than  50  knots 
s 

(Categories  B  and  C),  the  height  of  the  LMW  may  vary  widely  at  adjacent  grid  points 
and/or  stations.  This  is  a  common  characteristic  of  adjacent  wind  profiles  that  exhibit 
relatively  weak  maximums.  Thus,  it  is  not  surprising  that  the  rms  errors  for  Z(L) 
increase  with  decreasing  maximum  wind  speed. 

Conversely,  the  rms  errors  for  the  shears  below  and  above  the  LMW  decrease 
with  increasing  wind  speed,  but  have  approximately  the  same  percent  error  because 
the  values  of  vertical  wind  shear  are  generally  greater  for  the  higher  than  for  the  lower 
maximum  wind  speeds. 

The  total  rms  errors  for  individual  observation  times  for  the  initial-guess  and 
final  pass  for  each  of  the  four  LMW  parameters  is  graphically  represented  in  Figs.  6, 

7  and  8. 

9.  Objective  (computer)  versus  Subjective  (man)  LMW  Analysis  Comparison 

Although  the  areal-mean-error  method  of  verification  gives  useful  information 

about  the  quality  of  the  analysis  technique,  the  evaluation  is  not  complete  unless  a 

comparison  is  made  of  the  objective  computer  analysis  of  the  LMW  parameters  and 

the  subjective  analyses  performed  by  a  human  analyst.  Comparisons  for  each  of  the 

LMW  parameters  are  discussed  individually  below. 

Wind  Speed  at  the  LMW— Figure  9  illustrates  some  of  the  objective— subjective 

comparisons  of  the  jet  cores,  the  associated  maximum  isotachs,  and  the  75-knot 

isotachs.  These  and  the  other  objective  and  subjective  analyses  of  W  (L)  for  the  twelve 

s 

observation  times  reveal  that  there  are  only  minor  differences  in  the  analyses  over  the 
dense  data  areas,  while  more  differences  are  apparent  over  the  sparser  data  areas. 
These  differences  are  attributed  primarily  to  the  human  analyst  drawing  to  the  available 
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Fig.  6,  Rms  errors  of  wind  speed  at  the  LMW  for  three  wind-speed  categories  for  12  observation 
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Fig.  8,  Rms  errors  of  shears  below  (a)  and  above  (b)  the  LMW  for  three  wind-speed  categories  for  12 
observation  times. 


Fig.  9(a).  Objective  and  subjective  jet  stream  analysis  comparison,  1200Z, 
January  1,  1963  (knots). 
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Fig.  9(b).  Objective  and  subjective  jet  stream  analysis  comparison,  1200Z, 
January  5,  1963  (knots). 
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data,  resulting  in  a  tendency  for  the  jet  cores  to  be  positioned  close  to  the  observing 
station  when,  in  reality,  the  maximum  wind  speeds  may  be  further  from  the  station. 

A  subjective  analysis  could  not  be  performed  over  the  no-data  areas,  but  the 
nature  of  this  computer  analysis  technique  provides  for  an  analysis  over  the  entire  grid 
[e.g.,  the  initial  guess  of  the  LMW  wind  speed  uses  200-  and/or  300-mb  wind  speeds 
as  predictors  (see  Table  VI  and  the  Appendix)].  These  wind  speeds  are  derived  from  200- 
and  300-mb  height  analyses  by  using  the  gradient-wind  equation.  The  initial  guesses  to 
the  constant-pressure-surface  height  analyses  are  12-hr  d5mamic  prognoses  that  verify 
(terminate)  at  analysis  time.  Thus,  information  from  the  data  areas  (for  the  previous 
12  hrs)  is  incorporated  over  the  no-data  and  sparse-data  areas  of  the  LMW  analyses 
in  a  quantitative,  consistent,  and  objective  manner.  For  example,  the  ocean  areas  off 
the  east  coasts  of  continents  are  generally  no-data  areas,  but  the  dynamic  prognosis 
is  usually  based  on  the  initial  conditions  over  the  continent.  In  this  study,  the  300-  and 
200-mb  height  analyses  were  not  available  to  the  human  analyst.  Had  they  been  avail¬ 
able,  the  analyst  would  have  incorporated  the  information  in  a  subjective  manner  that 
probably  would  not  be  as  quantitatively  consistent  as  the  objective  procedure. 

Height  of  the  LMW— As  mentioned  previously,  the  analysis  of  Z(L)  is  only  meaning¬ 
ful  over  those  areas  where  W  (L)  >  75  knots.  Examination  of  the  objective  and  subjective 

s 

analyses  of  Z(L)  over  these  areas  reveals  that  they  are  in  general  agreement,  although 
the  objective  technique  has  a  tendency  to  indicate  cells  of  high  and  low  heights  near 
individual  stations. 

Figure  10  illustrates  the  objective  versus  the  subjective  Z(L)  charts  for  the  same 
observation  times  as  for  the  wind  speed  at  the  LMW.  In  some  cases  there  appear  to 
be  important  differences  in  the  analyzed  values  over  a  particular  area,  but  this,  pri¬ 
marily,  is  due  to  the  stations  that  were  withheld  from  the  objective  analysis  and  included 
in  the  subjective  analysis. 

Shears  Below  and  Above  the  LMW— The  objective  and  subjective  analyses  of 
vertical  shears  of  horizontal  wind  speed  below  and  above  the  LMW  exhibit  general 
patterns  that  are  similar.  The  detailed  patterns  appear  to  differ  in  some  cases,  especially 
near  the  jet  cores  where  the  subjective  analyst  tended  to  elongate  the  maximum  shear 
centers  while  the  objective  analysis  indicated  more  circular  maximums.  The  magnitude 
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Fig.  10(a).  Objective  and  subjective  height  of  LMW  analysis  comparison, 
1200Z,  January  1,  1963  (10^  ft). 
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Fig.  10(b).  Objective  and  subjective  height  of  LMW  analysis  comparison, 
1200Z,  January  5,  1963  (10^  ft). 
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of  the  shears  is  similar.  Figures  11  and  12  show  the  analyses  of  shear  below  and  above 
the  LMW  for  the  same  four  observation  times  as  for  the  wind-speed  and  the  wind-height 
analyses. 

10.  Consistency  Between  LMW  Parameters 

An  important  phase  of  the  evaluation  was  to  check  for  consistency  between  the 
LMW  parameters.  A  study  of  the  five  LMW  fields  for  each  of  the  observation  times  in 
the  January,  1963  data  sample  indicated  that  the  analyses,  in  general,  were  consistent 
with  one  another.  The  height  of  the  LMW  was  usually  at  a  higher  elevation  for  jets  at 
lower  latitudes  (subtropical  jets)  and  at  a  lower  elevation  for  mid-  and  higher -latitude 
jets  (polar  jets).  Where  two  jets  were  relatively  close  to  each  other  in  the  horizontal, 
higher  elevations  of  the  height  of  the  LMW  were  found  in  the  jet  to  the  right,  when  facing 
in  the  downstream  direction,  (i.e.,  the  direction  toward  which  the  wind  is  blowing). 

This  is  in  agreement  with  what  is  known  to  occur. 

The  maximum  values  of  vertical  shears  of  the  horizontal-wind  speed  frequently 
were  observed  to  be  in  the  regions  of  maximum  wind  speeds. 

Comparisons  between  Figs.  9  through  12,  and  between  the  five  parts  of  Fig.  13, 
will  verify  these  statements. 

The  exceptions  to  the  general  cases  occur  in  regions  where  the  vertical  wind 
profiles  are  blunt  in  the  region  of  maximum  wind  (i.e.,  where  there  is  really  no  clear 
level  of  maximum  wind  and  where  high  wind  speeds  exist  over  a  relatively  thick  layer 
of  8—15  thousand  feet  or  more  [4]).  Because  the  analysis  technique  specifies  a  level 
of  maximum  wind,  a  particular  level  is  indicated  somewhere  in  the  thick  layer  of  high 
wind  speeds,  although  the  level  selected  in  these  cases  is  not  meaningful  by  itself. 

The  shear  below  and  above  the  level  of  maximum  wind  is  averaged  for  layers 
that  are  10—12  thousand  feet  thick.  In  the  cases  of  thick  vertical  layers  of  high  wind 
speed,  either  the  shear  below,  above,  or  a  combination  of  both  may  be  computed  to  be 
a  relatively  low  value,  depending  on  where  in  the  layer  the  level  of  maximum  wind  is 
specified.  Solutions  to  the  problems  caused  by  thick  layers  of  maximum  wind  speeds 
are  offered  in  Section  V. 

The  following  interesting  items  should  be  noted  in  the  five  parts  of  Fig.  13; 

(a)  The  generated  observations  along  the  jet  cores  tend  to  elongate  the 


41 


Fig.  11(a).  Objective  and  subjective  shear  below  the  LMW  analysis  comparison, 
1200Z,  January  1,  1963  (knot  10~^  ft). 
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Fig.  11(b).  Objective  and  subjective  shear  below  the  LMW  analysis  comparison, 
1200Z,  January  5,  1963  (knot  10“^  ft). 
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Fig.  12(a).  Objective  and  subjective  shear  above  the  LMW  analysis  comparison 
1200Z,  January  1,  1963  (knot  10“3  ft). 
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Fig.  12(b).  Objective  and  subiective  shear  above  the  LMW  analysis  comparison, 
1200Z,  January  5,  1963  (knot  10“^  ft). 
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Fig.  13(a).  Objective  LMW  wind  speed,  initial  guess,  1200Z,  January  3,  1963 
(knots). 
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•  withheld  8tatio|ift 

■  Generated  observations  along  jet  core 

Fig.  13(b).  Objective  LMW  wind  speed,  final  pass,  1200Z,  January  3,  1963  (knots). 
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Fig.  13(c).  Objective  height  of  the  LMW  analysis,  1200Z,  January  3,  1963  (10^  ft). 
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Fig.  13(d).  Objective  shear  below  the  LMW  analysis,  1200Z,  January  3,  1963 
(knot  10"^  ft). 
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Fig.  13(e).  Objective  shear  above  the  LMW  analysis,  1200Z,  January  3  1963 
(knot  10-3  ft). 
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maximum  wind  isotachs.  Note  the  differences  between  the  initial-guess  analysis 
[Fig.  13(a)]  and  final -pass  analysis  [Fig.  13(b)],  particularly  over  the  areas  south 
of  Korea  and  north  of  the  Dominican  Republic.  The  jet  cores  were  drawn  taking 
into  account  the  wind  speed,  wind  direction,  and  generated  observations  along 
the  cores. 

(b)  The  subtropical  jet  (with  south- southwest  winds)  that  is  shown  in  the 
central  Pacific  [Fig.  13(c)]  meets  the  polar  jet  (with  west  winds)  in  the  vicinity 
of  50°  North  and  160°  West. 

(c)  The  jets  over  North  America  are  the  polar  jet  (elevation  28—32  thousand 
feet),  evident  from  central  Canada  southeastward  to  west  of  the  Great  Lakes,  and 
the  subtropical  jet  (elevation  36—40  thousand  feet)  which  appears  over  the  south 
central  United  States  and  runs  southeastward  into  the  southeastern  Gulf  of  Mexico, 
reappearing  south  of  Cuba  and  continuing  on  into  the  central  Atlantic. 

(d)  There  is  a  double  jet  in  the  western  Pacific,  with  the  northern  one  the 
polar  jet  (elevation  32—36  thousand  feet)  and  the  southern  one  the  subtropical  jet 
(elevation  36—40  thousand  feet). 

(e)  Over  most  areas,  the  maximum  shears  [Figs.  13(d)  and  (e)]  correspond 
closely  to  the  jet  cores. 
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SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 
The  following  conclusions  are  drawn: 

(a)  The  relatively-low  rms  errors  for  the  LMW  fields,  as  given  by  the  areal- 
mean-error  method,  indicate  the  analyses  are  of  high  quality. 

(b)  An  advantage  of  the  objective  technique  over  the  subjective  analysis  is 
that  it  provides  analyses  over  areas  of  the  hemisphere,  including  no-data  and 
sparse-data  areas.  Information  from  data  areas  was  incorporated  into  the  LMW 
analyses  over  the  sparse-  and  no-data  areas  in  a  quantitative,  objective,  and 
consistent  manner  by  using  LMW  screening  regression  equations  for  seven  jet- 
stream  categories  in  conjunction  with  constant-pressure-surface  fields. 

(c)  Over  the  dense  data  areas,  the  objective  analyses  compare  favorably 
with  the  subjective  analyses. 

(d)  The  objective  analyses  for  the  LMW  parameters  are  consistent  with  one 
another. 

(e)  The  generated  jet  core  “observations”  tend  to  elongate  the  maximum 
isotachs. 

(f)  An  apparent  deficiency  in  the  model  is  that  vertical  wind  profiles  exhibit¬ 
ing  a  thick  layer  of  high  wind  speeds  are  not  specified  in  the  analysis  of  the  LMW 
parameters.  Because  the  analysis  technique  requires  a  level  of  maximum  wind, 

a  particular  level  is  selected  somewhere  in  the  thick  layer,  although  the  level 
selected  is  not  meaningful  for  these  cases.  The  analyses  and  verifications  of 
Z(L),  S^,  and  are  sensitive  to  where  the  level  of  maximum  wind  is  specified 
in  the  layer. 

Based  upon  these  conclusions,  the  following  is  recommended: 

(a)  A  layer  of  maximum  wind  should  be  defined  similar  to  that  discussed  by 
Reiter  [4],  with  one  of  the  exceptions  being  that  the  analysis  of  wind  speed  should 
be  for  the  peak  wind  speed  in  the  layer,  rather  than  for  the  mean  maximum  in  the 
layer  that  Reiter  defines  as  90%  of  the  peak  wind.  The  reason  for  analyzing  for 
100%  of  the  peak  wind  is  that  the  LMW  analysis  is  applied  to  clear-air  turbulence 
analysis  and  forecasting. 
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(b)  The  mean  altitude  of  the  layer  of  maximum  wind,  rather  than  the  height 
of  the  level  of  maximum  wind,  should  be  analyzed.  The  thickness  of  the  layer  of 
maximum  wind  should  be  analyzed  as  an  additional  field  that  will  supply  more 
detailed  information  about  the  three-dimensional  structure  of  the  wind  field  in 
the  vicinity  of  the  jet  stream. 

(c)  To  complete  the  more  realistic  three-dimensional  picture,  the  vector 
shear  of  the  wind  below  and  above  the  layer  of  maximum  wind  should  be  analyzed 
instead  of  the  scalar  shear  below  and  above  the  level  of  maximum  wind. 

To  accomplish  the  work  suggested  by  these  recommendations,  it  will  be  necessary 
to;  (1)  redefine  some  of  the  jet  stream  categories  by  considering  the  thickness  of  the 
layer  of  maximum  wind  as  an  additional  criterion  in  the  category  definitions,  (2)  process 
the  station  data  according  to  the  LMW  fields  listed  in  the  recommendations,  (3)  derive 
new  LMW  regression  equations,  based  upon  the  new  category  definitions,  for  specifying 
initial  guesses  to  the  LMW  parameters,  and  (4)  add  the  thickness  of  the  layer  of  maxi¬ 
mum  wind  to  the  LMW  analysis  program. 
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EQUATIONS  FOR  LMW- 
PARAMETER  FIRST-GUESS  FIELDS 


APPENDIX 


EQUATIONS  FOR  LMW- PARAMETER  FIRST-GUESS  FIELDS 

1.  Category  1— No  Jet 

Conditions: 

W  (3)  and  W  (2)  <  50  knots 
s  s 

Z(L)  =  2908.9  +  49.977[W  (1)]  -  37.794[W  (7)]  +  23.941[T(7)] 

s  s 

W  (L)  =  14.190  +  0.36354[W  (2)]  +  0.36641[W  (1)]  +  0.30954[W  (5)] 

s  s  s  s 

S  =  0.013107  -  0.14963(S  )  +  0.56996[W  (1)]  +  0.067427[W  (L)] 

3.  OSS 

S  =  2.49333  -  0.064107[W  (L)]  +  0.030697[W  (2)]  +  0.012455[W  (5)] 

o  s  s  s 

+  0.00022536[Z(L)] 

2.  Category  2— Weak  to  Moderate  Subtropical  Jet 

C  onditions : 

(a)  W  (3)  or  W  (2)  =  50—99  knots  with  latitude  <  30°N,  or 

s  s 

(b)  W  (2)  -W  (3)  >  10  knots  with  latitude  45°N 

s  s 

Z(L)  =  4048.8  +  12.306[W  (1)]  -  8.0158[W  (3)]  +  12.709[T(5)] 

s  s 

W  (L)  =  13.504  +  0.94295[W  (2)] 

s  s 

S  =  0.80542  -  0.12514[W  (L)J+  0.04950[W  (1)]+  0.05469[W  (2)] 

3,  S  S  S 

S  =  13.439  -  0.0037641[Z(L)]  +  0.13056[W  (L)]  -  0.10909[W  (3)] 

D  S  S 

3.  Category  3— Strong  Subtropical  Jet 

Conditions: 

(a)  W  (3)  or  W  (2)  S;  lOO  knots  with  latitude  <  30°N,  or 

S  S 

(b)  W  (2)  -W  (3)  >  15  knots  with  latitude  45°N 

S  s 

Z(L)  =  3822.1  -  4.8446[W  (3)]  +  3.0291[W  (1)]  +  2.7354[W  (2)] 

s  s  s 

W  (L)  =  27.774  +  0.84052[W  (2)] 

s  s 

S^  =  10.846  -  0.15284[Wg(L)]  +  0.087687[Wg(2)]  -  0.0026408[Z(L)] 

+  0.054073[W  (1)] 
s 
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S,  =  17.851  -  0.0046131[Z(L)]  -  0.098653[W  (3)]  +  0.10677[W  (L)] 

u  S  S 

4.  Category  4— Weak  to  Moderate  Subtropical  Polar  Jet 

Conditions: 

W  (3)  or  W  (2)  =  50-99  knots,  with  IW  (2)  -  W  (3)1  10  knots  and 

s  s  s  s 

latitude  =  30— 45°N 

Z(L)  =  3604.9  +  11.855[W  (1)]  +  16.843[T(7)]  -  7.4020[W  (5)]  -  34.071[W  (3)] 

s  s  s 

+  32.035[W  (2)] 
s 

W  (L)  =  13.623  +  0.97940[W  (2)] 

s  s 

S  =  0.35951  -  0.13492[W  (L)]  +  0.089316[W  (3)]  +  0.038302[W  (1)] 

3.  S  S  S 

S  =  8.5199  +  0.15191[W  (L)]  -  0.  002418[Z(L)]  -  0.10463[W  (3)]  -  0.05478[W  (5)] 

u  S  S  S 

5.  Category  5— Strong  Subtropical  Polar  Jet 

Conditions: 

W  (3)  or  W  (2)  s  100  knots  with  IW  (2)-W  (3)1  15  knots 

s  s  s  s 

and  latitude  30°— 45°N 

Z(L)  =  3727.3  +  11.336[T(7)]  +  5.3228[W  (1)]  -  22.056[W  (3)]  +  18.295[W  (2)] 

s  s  s 

W  (L)  =  21.365  +  0.48576[W  (3)]  +  0.45093[W  (2)] 

s  s  s 

S  =  9.0880  -  0.0032362[Z(L)]  -  0.35472(8  )  +  0.04686[W  (1)]  -  0.080768[W  (L)] 

D  S  S 

+  0.066046[W  (2)] 
s 

S  =  12.340  -  0.0038816[Z(L)]  +  0.14103[W  (L)]  -  0.11350[W  (3)] 

o  s  s 

6.  Category  6— Weak  to  Moderate  Polar  Jet 

Conditions: 

(a)  W  (3)  or  W  (2)  =  50—99  knots  with  latitude  >  45°N,  or 

s  s 

(b)  W^(3)  -W^(2)  >10  knots  with  latitude  >  30°N 

Z(L)  =  3361.9  -  10.355[W  (5)]  +  16.909[W  (2)]  -  10.463[W  (3)] 

s  s  s 

W  (L)  =  17.319  +  0.64497[W  (3)]  +  0.30167[W  (2)] 
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S  =  2.2800  -  0.090459(W  (L)]  +  0.054284(W  (2)]  -  0.00082532[Z(L)] 

3,  S  S 

+  0.033011[W  (1)] 

S 

S  =  4.8737  +  0.10729[W  (L)]  -  0.065994[W  (5)]  -  0.0015093[Z(L)] 

o  s  s 

-  0.038469(W  (3)] 
s 

7.  Category  7— Strong  Polar  Jet 

Conditions: 

(a)  W  (3)  or  W  (2)  ^  100  knots  with  latitude  >  45°N,  or 

s  s 

(b)  W  (3)  -W  (2)  >  15  knots  with  latitude  ^  30°N 

s  s 

Z(L)  =  3794.5  +  20.950[T(5)]  +  5.9493[W  (2)]  -  5.6606[W  (3)] 

s  s 

W  (L)  =  31.007  +  0.64466[W  (3)]  +  0.23493[W  (2)] 

s  s  s 

S  =  8.6349  -  0.16002[W  (L)]  +  0.090277[W  (2)]  -  0.0028288[Z(L)] 

3,  S  S 

+  0.060576[W  (3)] 
s 

S  =  11.040  -  0.0035587(Z(L)]  +  0.083069[W  (L)]  -  0.083452(W  (5)] 

o  s  s 
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SECTION  I 


INTRODUCTION 

Two  objective  techniques  for  predicting  level  of  maximum  wind  (LMW)  parameters 
were  developed.  The  first,  the  LMW  Modeling  Prediction  Technique,  consists  simply  of 
applying  the  diagnostic  equations  that  specify  LMW  parameters  (described  in  PART  I) 
to  prognoses  of  constant-pressure-surface  fields.  This  results  in  prognoses  of  the 
following  parameters: 

(a)  wind  speed  at  the  LMW  [W^(L)] 

(b) '  wind  direction  at  the  LMW  [W^(L)] 

(c)  height  of  the  LMW  [Z(L)] 

(d)  shears  above  and  below  the  LMW  [S  and  S,  ] 

'  '  a  b 

(In  PART  I,  the  equations  were  used  with  observed  fields  to  generate  initial  guesses  for 
the  LMW  analyses.) 

An  additional  feature  of  this  technique  is  that  it  locates  the  forecast  position  of 
the  jet  core,  generates  “observations”  along  the  core  using  horizontal  profiles  of  jet- 
stream  models  [1],  and  adjusts  the  wind-speed  grid  with  the  successive  approximation 
analysis  technique  (SAT)  and  the  generated  observations. 

This  technique  is  simply  a  by-product  of  the  LMW  analyses  (PART  I),  whereas  the 
following  technique  is  designed  specifically  for  predictions. 

The  second  technique,  the  LMW  Predictor  Network  Technique  (PNT),  is  more 
complex  and  provides  12-  and  24-hr  prognoses  of  LMW  parameters  using  a  moving- 
coordinate  prediction  system.  For  PNT,  regression  equations  were  derived,  for  12-  and 
24-hr  prediction  of  the  parameters  or  their  changes,  for  grid  points  categorized  accord¬ 
ing  to  the  initial-state  (t^)  position  of  the  jet  stream  [as  indicated  by  the  analysis  of 
W^(L)  and  W^(L)]. 

The  possible  predictors  available  for  selection  in  the  screening  procedure  were 
the  five  LMW  parameters  at  t^,  their  previous  12-hr  changes,  the  wind  speed  at  500, 

300,  200,  and  100  mb  (at  t^),  the  wind-speed  shears  for  the  300-500-,  200-300-,  and 
100— 200-mb  levels,  and  their  previous  12-hr  changes.  These  predictors  were  selected 
for  screening  by  applying  a  grid  network  oriented  with  the  jet-stream  flow  to  the  area 
about  the  predictand  point. 
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SECTION  II 


DATA  PROCESSING  FOR  THE  PREDICTOR  NETWORK  TECHNIQUE 

Hemispheric  grid-point  LMW  analyses  (see  lART  I  of  this  report)  were  generated 
from  data  (2  observations  per  day)  provided  by  the  National  Meteorological  Center  (NMC) 
for  20  days,  of  both  the  months  of  December  1963  and  January  1964.  Wind  speed 
analyses  for  the  four  constant-pressure  levels  used  in  the  study  were  available  on  mag¬ 
netic  tape  from  Global  Weather  Central  (GWC).  The  developmental  sample  consisted 
of  the  first  fifteen  days  of  data  of  each  month.  The  remaining  five  days  were  used  as 
independent  data  for  testing. 

After  all  basic  analysis  grids  were  obtained,  the  predictand  points  were  selected 
on  the  basis  of  objective  criteria  described  in  Section  III.  The  predictand  points  were 
selected  only  within  the  a,rea  shown  in  Fig.  1,  as  the  main  purpose  of  the  feasibility 
study  was  to  test  a  prediction  technique  that  would  be  applicable  in  the  area  about  North 
America.  The  grid  network  was  then  applied  to  each  predictand  point,  and  oriented 
according  to  the  W^(L)  at  the  predictand  point;  the  predictor  information  from  each  field 
was  extracted. 
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Fig.  1.  Area  used  for  North  American  jet-stream  prediction  study.  Latitude 
33°N  is  indicated  as  the  dividing  line  between  polar  and  subtropical  jet  stream 
categories. 


SECTION  III 


THE  LMW  PREDICTOR  NETWORK  TECHNIQUE 

1.  Categorization  of  Predictand  Points 

The  development  of  this  technique  to  predict  the  LMW  parameters  required  the 
careful  utilization  of  all  the  initial-state  information  that  describes  the  parameter 
fields.  A  procedure  in  which  prediction  equations  are  derived  for  application  at  all 
grid  points,  such  as  one  used  in  the  prediction  of  100-mb  height  change  [3],  did  not 
seem  a  fruitful  approach  here  because  of  the  characteristics  of  the  predictands.  Per¬ 
haps  the  most  important  feature  of  any  LMW  prediction  technique  is  the  forecasting  of 
the  location  of  the  jet-stream  core  and  the  corresponding  wind-speed  variations,  or  jet- 
maxima,  located  along  the  core.  The  variations  in  the  fields  of  vertical  wind  shear  above 
and  below  the  LMW  and  the  height  of  the  LMW  are  closely  related  to  the  position  of  the  jet- 
stream  core.  It  followed  from  this  that  prediction  equations  should  be  developed  for 
application  to  grid  points  which  are  first  categorized  relative  to  the  initial- state  position 
of  the  core.  (The  procedure  for  determining  the  core  and  selecting  those  consecutive 
Numerical  Weather  Prediction  (NWP)  grid  points  that  most  closely  depict  the  core  is 
given  in  PART  1.) 

The  procedure  for  categorizing  the  NWP  grid  points,  relative  to  the  analyzed  position 
of  the  core  at  time  t^,  is  illustrated  in  Fig.  2.  In  this  figure,  a  solid  line  connects  those 
consecutive  grid  points  that  most  closely  approximate  the  analyzed  core.  The  arrow 
indicates  the  predominant  wind  direction  along  the  core.  Five  categories  are  established. 
The  grid  point  labeled  “A”  and  circled  indicates  the  approximate  location  of  the  jet-stream 
maximum.  The  wind  speed  at  this  grid  point  is  greater  than  that  at  the  adjacent  core  points 
upstream  (the  direction  from  which  the  wind  is  blowing)  or  downstream.  The  grid  point 
labeled  “B”  is  two  core  points  upstream  and  the  one  labeled  “C”  is  two  core  points  down¬ 
stream.  Grid  point  “D”  is  the  NWP  grid  point  closest  to  two  grid  intervals  to  the  left  of  the 
core  (using  the  convention  of  looking  downstream)  and  grid  point  “E”  is  closest  to  two  grid 
intervals  to  the  right. 

The  above  approach  establishes  five  categories  of  grid  points  relative  to  a  given 
jet-stream  system.  However,  there  are  at  least  two  basic  types  of  jet-stream  systems— 
polar  and  subtropical— and  the  characteristics  and  behavior  of  these  two  systems  differ 
considerably.  For  example,  the  subtropical  jet  often  persists  in  the  same  general  location 
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Fig.  2.  Procedure  for  categorizing  NWP  grid  points  from  the 
initial-state  position  of  the  jet  stream  core.  Solid  line  connects  the 
core,  and  the  categorized  grid  points  are  marked  (see  text  for 
explanation) . 

with  predominant  west  to  east  flow  that  varies  only  moderately  as  jet-maxima  move 
through  the  system.  The  winds  in  the  polar  jet,  on  the  other  hand,  often  have  a  significant 
north-south  component  of  flow,  develop  and  dissipate  rapidly,  and  experience  rapid 
changes  in  location  over  short  periods  of  time.  Because  of  the  important  behavioral 
differences  in  the  two  jet-stream  types,  we  developed  separate  prediction  relationships 
for  the  LMW  parameters  based  on  whether  the  jet-stream  core  was  subtropical  or  polar 
(see  Table  I  for  the  grid  point  categories). 

It  is  obviously  difficult  to  develop  objective  criteria  for  classifying  the  jet  core 
as  polar  or  subtropical.  For  this  limited  feasibility  study,  a  jet  core  was  considered 
polar  if  grid  point  A  (wind  speed  maximum)  of  the  classification  model  shown  in  Fig.  2 
was  located  north  of  33°N,  and  subtropical  if  located  south  of  this  latitude.  It  is  recog¬ 
nized  that  a  separation  of  the  two  jet  types  based  only  on  latitude  is  rather  crude.  How¬ 
ever,  we  felt  that  a  more  sophisticated  scheme  was  not  justified  in  the  first  attempt  to 
examine  this  forecasting  problem. 
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TABLE  1 

GRID  POINT  CATEGORIZATION 


Polar 

category 

Subtropical 

category 

Definition 

- 1 

lA 

1  '  - 

1  2A 

Wind  speed  maximum 

IB 

2B 

Second  core  point  upstream* 

1C 

2C 

Second  core  point  downstream* 

ID 

2D 

Two  grid  intervals  to  the  left  of  coref 

IE 

2E 

Two  grid  intervals  to  the  right  of  coret 

^Upstream  is  direction  from  which  wind  is  blowing.  Downstream  is 
the  direction  toward  which  wind  is  blowing, 

t Convention  used  is  that  of  looking  downstream. 

2.  Predictor  Network  Grid 

The  predictor  network  grid  used  to  select  predictors  at  NWP  grid  points  is  shown 
in  Fig.  3.  The  grid  interval  is  one  NWP  grid  interval  (381  km  at  60”N).  NWP  grid- 
point  analyses  of  all  predictand  and  predictor  fields  were  contained  on  magnetic  tape. 

The  predictand  points  are  selected  by  the  objective  procedure  described  in  Section  1. 

The  predictor  network  grid  is  then  applied  to  all  predictand  points.  The  values  of  the 
predictors  are  determined  for  all  grid  points  of  the  network  by  a  computer  program. 

The  grid  is  placed  on  each  predictor  field  such  that  the  predictand  grid  point  is  located 
at  K-5,  L=3,  and  the  L  axis  of  the  grid  is  oriented  parallel  to  the  LMW  wind  direction  at 
location  (5,  3)  (see  Fig.  3).  Other  orientations  of  the  network  grid,  such  as  the  L  axis 
always  being  east— west  or  always  parallel  to  the  basic  NWP  grid,  were  possible  but 
we  decided  that  the  sampling  of  the  predictor  fields  at  locations  relative  to  the  pre¬ 
dictand  point  would  be  considerably  more  meaningful  if  the  network  was  oriented 
relative  to  the  flow  of  the  jet-stream  system. 

A  7x5  (rather  than  a  square)  array  was  used  so  that  rapidly  moving  features  well 
upstream,  such  as  jet -maxima  or  short-wave  troughs,  would  come  within  the  a.rea  defined 
by  the  predictor  network  grid. 
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Fig.  3.  Predictor  network  grid.  Predictand  point  is  at  K=5,  L=3. 
[The  predictor  points  are  indicated  by  appropriate  symbol  (see  Section 
4. for  meaning)]. 


3 .  Predictands 

The  predictands  (the  parameters  being  forecast)  considered  in  this  study  were 
the  12-  and  24-hr  changes,  and  the  actual  values,  of  the  five  LMW  parameters  [wind 
direction,  W  (L);  wind  speed,  W  (L);  height,  Z(L);  and  the  shear  above  and  below  the 

Q  S 

LMW,  S  and  S  ].  This  gives  a  total  of  twenty  predictands  for  each  of  the  ten  categories 

Si  D 

described  in  Section  1.  The  predictands  are  shown  in  Table  II.  Obviously,  selectivity 
was  required  in  the  feasibility  study  experiments.  Greatest  emphasis  was  placed  on 
predicting  the  12-  or  24-hr  changes  of  the  predictands,  not  on  predicting  their  actual 
value  12  or  24  hours  in  the  future. 

4.  Predictors 

The  7x5  network  grid  is  applied  to  20  fields  of  possible  predictors  that  have  been 
analyzed  on  a  NWP  hemispheric  grid.  If  every  point  on  the  grid  array  for  each  field 
were  utilized,  there  would  be  a  total  of  700  possible  predictors. 

Because  the  screening-regression  technique  of  Miller  [2],  as  programmed  for  the 
IBM  7094  computer,  has  an  upper  limit  of  180  possible  predictors  and  predictands,  it 
was  necessary  to  select  subjectively  from  the  original  700  predictors,  a  set  of  160  or 
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TABLE  II 

POSSIBLE  PREDICTANDS 


Predictand* 

Forecast 
interval  (hr) 

Predictand 
symbol  t 

12 

W^(L) 

24 

*d<^>24 

12 

^*d<‘-'l2 

AW,(L) 

a 

24 

^"'d<‘->24 

W  (L) 

12 

s 

24 

*s<‘->24 

12 

™s<'‘>12 

AW  (L) 

s 

24 

^*s<^'24 

12 

Z(L) 

24 

AZ(L) 

12 

AZ(L)^2 

24 

AZ(L)2^ 

12 

®bl2 

s. 

b 

24 

®b24 

AS^ 

12 

^"b  12 

b 

24 

^^b24 

12 

®al2 

S 

a 

24 

\  24 

AS 

12 

“a  12 

a 

24 

^"a24 

*DeIta  (A)  indicates  ‘‘change  in.” 

tThe  hat  (- )  is  used  to  identify  a  predictand. 

A  parameter  symbol  without  the  hat  would  indicate 
a  predictor. 
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less  for  screening.  This  restriction  is  not  very  severe  when  one  considers  the  high 
correlation  that  exists  between  values  of  predictors  from  the  same  field  at  adjacent 
grid  points. 

The  number  and  location  in  the  network  grid  of  screened  predictors  from  each 
field  can  be  seen  by  examining  Fig.  3.  and  Table  III.  The  three  types  of  distributions  of 
predictors  available  for  the  screening-regression  technique  are: 

(a)  Selection  Type  A— all  grid  points  in  Fig.  3  enclosed  by  a  circle, 
square  or  triangle. 

(b)  Selection  Type  B— all  grid  points  in  Fig.  3  enclosed  by  a  circle  or 
square . 

(c)  Selection  Type  C— all  grid  points  in  Fig.  3  enclosed  by  a  square. 

Table  III  lists  each  predictor  field,  distribution  selection  type,  and  the  number  of 
predictors  screened. 
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TABLE  III 

PREDICTORS  AVAILABLE  FOR  SCREENING 


Predictor 

Description 

Number 

available 

Selection 

type 

W^(L) 

LMW  wind  direction 

6 

c 

VV  (L) 
s 

LMW  wind  speed 

18 

A 

Z(L) 

LMW  height 

6 

C 

"b 

Shear  below  LMW 

12 

B 

s 

a 

Shear  above  LMW 

12 

B 

W  (5)* 
s 

500-mb  wind  speed 

6 

C 

W^(3) 

300-mb  wind  speed 

6 

C 

Ws(2) 

200-mb  wind  speed 

6 

C 

W^(l) 

100-mb  wind  speed 

6 

C 

AWd(L)_i2 

12-hr  LMW  wind  direction  change 

1 

Predictand 
point  only 

AWs(L)_i2 

12-hr  LMW  wind  speed  change 

12 

B 

^Z(L)_12 

12-hr  LMW  height  change 

6 

C 

^^-12 

12-hr  shear  below  LMW  change 

12 

B 

^^a-12 

12-hr  shear  above  LMW  change 

12 

B 

S  (3  -5) 

300— 500-mb  wind  speed  difference 

6 

C 

S(2  -3) 

200— 300-mb  wind  speed  difference 

6 

C 

S(l-2) 

100— 200-mb  wind  speed  difference 

6 

c 

AS  (3-5)_^2 

12-hr  change  of  S  (3  -5) 

6 

c 

AS  (2  -3)_^2 

12-hr  change  of  S(2-3) 

6 

c 

12-hr  change  of  S(l-2) 

6 

c 

^Parenthetic  numbers  indicate  level  (in  hundreds  of  mb). 


SECTION  IV 

RESULTS  USING  THE  PREDICTOR  NETWORK  TECHNIQUE 

A  series  of  initial  experiments  were  conducted  to  reduce  the  number  of  pre- 
dictands  and  predictors  used  in  the  derivation  of  prediction  equations  for  the  ten  grid 
point  categories.  This  experimentation  was  undertaken  in  an  attempt  to  conserve  com¬ 
puter  time,  and  because  it  is  obvious  that  both  the  predictands  and  the  predictor  fields 
contain  redundant  information.  Basically,  two  sets  of  experiments  were  performed. 

In  the  first  set  (reduction  of  predictands),  the  same  set  of  predictors  used  to  obtain 
regression  equations  for  the  future  12-  and  24-hr  values  of  the  LMW  parameters  were 
also  used  to  obtain  equations  for  the  12-  and  24-hr  changes  of  these  parameters.  In 
the  second  set  of  experiments  (reduction  of  predictors)  various  sets  of  predictors 
were  used  to  develop  regression  equations  to  forecast  the  same  predictand.  Both  sets 
of  experiments  were  undertaken  using  mainly  grid-point  categories  ID,  IE,  and  2C, 
which  contain  336,  343  and  104  dependent  cases,  respectively.  These  categories  were 
tested  most  extensively  because  their  associated  data  samples  were  first  available. 

Since  the  results  obtained  from  these  categories  and  limited  testing  with  other  cate¬ 
gories  was  similar,  extensive  preliminary  testing  with  all  categories  was  not  considered 
necessary.  The  detailed  results  of  all  the  initial  experiments  are  not  given  here. 

Rather,  a  summary  of  the  results  that  led  to  the  final  selection  of  predictands  and 
predictors  is  given. 

Table  IV  gives,  for  the  first  set  of  experiments,  the  total  percent  reduction  of 
variance  (%  red.)  of  the  LMW  predictands  for  grid-point  Category  IE.  The  LMW  wind 
direction  is  not  included  in  the  experiment.  It  was  found  that  the  skill  in  predicting 
wind  direction  was  very  low  and,  as  good  prognoses  of  wind  direction  are  presently 
available  from  NWP,  we  decided  to  discontinue  experiments  containing  wind  direction 
as  a  predictand. 

The  identical  predictors  were  screened  for  each  predictand  pair  (shown  in 
Table  IV)  (12-hr  LMW  wind  speed  and  12-hr  change  of  LMW  wind  speed  are  a  pre¬ 
dictand  pair).  The  predictive  skill,  as  evident  in  the  %  red.,  is  much  greater  in  fore¬ 
casting  12-  and  24-hr  changes  of  LMW  height,  shear  below,  and  shear  above  than 
in  forecasting  their  actual  future  values.  However,  in  the  case  of  the  LMW  wind  speed, 
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TABLE  IV 

REDUCTION  OF  PREDICTANDS 

[Comparison  of  percent  reductions  of  LMW  predictands 
using  same  predictor  set  (Category  IE)] 


%  red. 

Predictand 

12- hr 
prediction 

24-hr 

prediction 

W  (L) 
s 

57 

52 

AW  (L) 
s 

33 

36 

Z(L) 

34 

27 

AZ(L) 

43 

51 

19 

10 

48 

38 

s 

a 

29 

24 

AS 

a 

45 

46 

a  larger  %  red.  is  obtained  using  the  actual  value,  rather  than  the  change,  for  both  12- 
and  24-hr  prediction  intervals.  The  summary  results  shown  in  Table  IV  are  similar 
to  those  obtained  for  Categories  ID  and  2C,  except  that  all  values  of  %  red.  were  higher 
for  Category  2C  due  to  the  greater  persistence  of  the  subtropical  jet  stream. 

It  is  certainly  not  surprising  that  the  predictors  are  better  related  (higher  linear 
correlations)  to  the  change  of  height  or  wind  shears  rather  than  to  the  actual  values. 
For  example,  a  LMW  height  value  of  35000  ft  could  correspond, to  a  wide  variety  of 
meteorological  conditions,  while  an  LMW  height  change  of  -3500  ft  would  frequently 
be  associated  with  an  approaching  trough  or  wind  speed  maximum.  In  the  case  of  wind 
speed,  however,  the  actual  value  clearly  defines  the  jet  stream  and  associated  core 
maximum  and  is,  therefore,  well  correlated  to  potential  predictors.  On  the  basis  of 
the  initial  testing  of  the  three  categories,  we  decided  to  use  the  actu al  value  of  LMW 
wind  speed  and  the  changes  of  the  LMW  height  and  wind  shear  fields  as  predictands  to 
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develop  12-  and  24-hr  prediction  equations  for  the  ten  grid-point  categories. 

An  example  of  the  results  from  the  second  set  of  experiments  is  shown  in  Table  V 
which  gives,  for  Category  2C,  the  total  percent  reduction  of  variance  of  the  12-  and 
24-hr  LiMW  wind  speed,  and  changes,  when  (1)  all  predictors  are  screened  and  (2)  the 
fields  containing  wind-speed  differences  between  constant  pressure  surfaces,  and  the 
12-hr  changes  of  these  differences,  are  excluded.  Note  that  the  change  in  %  red.  is  very 
small  when  the  predictor  fields  S(3  -5),  S(2-3),  and  S(l-2),  and  their  12-hr  changes,  are 
eliminated  from  consideration  in  obtaining  the  regression  equations.  On  the  basis  of 
the  experiments  summarized  in  Table  V,  and  many  similar  experiments,  we  decided  to 
exclude  the  prediction  fields  containing  the  actual  values  and  the  previous  12-hr  changes 
of  the  wind- speed  differences  in  the  development  of  the  regression  equations  for  the  ten 
grid -point  categories. 

Various  other  initial  prediction  experiments  were  attempted,  with  different  com¬ 
binations  of  predictor  fields  being  included  for  screening.  The  general  conclusion  was 
that  the  actual  values  and  12-hr  changes  of  the  LMW  predictors,  and  wind  speed  at 
constant  pressure  surfaces,  all  contain  significant  predictive  information,  the  relative 
importance  of  which  depends  on  the  category,  forecast  interval  and  LMW  predictand. 
These  predictors  were  used  to  develop  the  prediction  equations  for  all  categories. 

TABLE  V 

REDUCTION  OF  PREDICTORS 

[Comparison  of  percent  reduction  of  Wg(L)  for 
different  predictor  sets  (Category  2C)] 


1 

Predictand 

%  red.  using 

All 

predictors 

All  except  S(3  -5),  S(2  -3) 
S(l-2),  and  their  changes 

Ws(L)i2 

76 

73 

AW^(L),2 

52 

47 

Ws(L)24 

59 

59 

AWg(L)24 

46 

45 
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5.  Polar  Jet  Stream— Dependent  Data  Results 

The  dependent  sample  size  for  each  of  the  five  polar  jet-stream  categories  was 
approximately  350  cases.  The  dependent  sample  was  obtained  from  54  observation 
times  within  the  two  intervals  12Z  Dec.  1 — 12Z  Dec.  14,  1963  and  OOZ  Jan.  2  — 

OOZ  Jan.  15,  1964,  Table  VI  lists  the  predictors  selected  for  each  of  the  eight  pre- 
dictands  for  each  grid-point  category,  and  the  total  percent  reduction  of  variance  (%  red.) 
of  the  predictand  for  the  predictors  selected.  In  examining  the  lists  of  variables  selected, 
the  most  obvious  repetitive  feature  is  that  the  first  predictor  selected  for  both  12-  and 
24-hr  prognoses  of  the  four  LMW  parameters  was,  in  most  instances,  the  initial-state 
value  at  the  predictand  point  (5,3)  of  the  LMW  parameter  being  predicted.  For  example, 
in  all  five  categories,  for  12-  and  24-hr  predictions  of  height  change,  the  first  pre¬ 
dictor  selected  was  Z(L.)  (5,3). 

A  second  feature  to  be  noted  is  that  very  few  12-hr  change  predictors  were  selected. 
This  indicates  that  the  12-hr  changes  prior  to  t^  (initial  state)  are  not  highly  correlated 
(linearly)  to  the  12-  and  24-hr  LMW  predictands.  This  again  reflects  the  fact  that  the 
polar  jet  streams  are  rapidly-changing  (non -linearly)  systems.  It  will  be  seen  that  the 
12-hr  change  predictors  contain  more  useful  information  when  forecasting  the  subtropical 
jet-stream  LMW  predictands. 

Two  additional  general  comments  can  be  made.  First,  aside  from  the  predictand 
point  (5,  3),  the  selected  predictors  are  most  frequently  located  at  a  point  upstream  from 
the  predictand  point.  (The  parenthetic  numbers  separated  by  a.  comma,  given  with  each 
selected  predictor,  indicate  the  location  on  the  predictor  grid.  When  the  first  of  the  two 
numbers  is  between  1  and  4,  an  upstream  location  is  indicated;  the  numbers  6  and  7 
indicate  a  downstream  location).  This  is  true  for  all  ca.tegories,  and  is  certainly 
reasonable  because  it  would  be  expected  that  most  of  the  useful  predictor  information 
would  be  found  upstream  (particularly  in  the  case  of  rapidly  changing  jet- stream  systems). 
Also  note  that  the  type  of  predictor  most  frequently  selected  is  the  wind  speed  either  at 
the  LMW  or  a,  constant  pressure  surface.  This  result  is  also  quite  reaonsable  because, 
of  the  predictor  fields  considered  in  the  feasibility  study,  the  wind  speed  most  clearly 
defines  the  jet-stream  systems,  and  therefore  contains  the  most  useful  prediction  infor¬ 
mation. 
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TABLE  VI 

SELECTED  PREDICTORS  FOR  POLAR  JET-STREAM  CATEGORIES 
(Results  on  dependent  data) 

(a)  Category  LA  (350  cases) 


Predictand 

Predictors  in  order  of  selection 

% 

red. 

Std. 

dev. 

rms 

error 

Units 

W^(2)(3.3).  W^(L)(4.2).  AS^(7. 1) 

40.7 

33.2 

25.5 

knots 

AZ(L)^2 

Z(L)(5,3),  AZ(L)(5,3) 

27.0 

5430 

4641 

ft 

^^bl2 

S^(5.3).  Wg(L)(4.2) 

52.9 

2.7 

1.8 

knots  10  ^  ft  ^ 

^®al2 

S  (5.3).  W  (L)(2.2) 

cL.  b 

53.2 

2.8 

1.9 

-3  -1 

knots  10  ft 

W^(L)(3.3).  W^(2)(3.1).  S^^(3.3) 

34.9 

32.4 

26.1 

knots 

AZ(L)24 

Z(L)(5.3).  Z(L)(1.3).  W  (1)(5.3) 

O 

30.8 

7043 

5860 

ft 

^^^b24 

Sb(5.3).  W^(2)(3.1) 

64.1 

2.6 

1.6 

knots  10  ^ft  ^ 

^^^24 

S^(5.3).  W^(2)(3.1) 

56.9 

2.9 

1.9 

knots  10  ^ft  ^ 

(b)  Category  IB  (355  cases) 


Predictand 

Predictors  in  order  of  selection 

% 

red. 

Std. 

dev. 

rms 

error 

Units 

W  (L)(3.3).  W  (5)(5.3).  W  (L)(6.  2). 

S  S  5 

W  (1)(5.3) 

D 

48.7 

33.3 

23.8 

knots 

AZ(L)^2 

Z(L)(5.3).  Z(L)(3.3).  W  (2)(3. 1). 

5 

Wg(5)(5.3).  Wg(L)(5.5) 

37.5 

5858 

4631 

ft 

S^j(5.3).  Wg(L)(6.2) 

59.2 

2.6 

1.6 

knots  10~^  ft  ^ 

S^(5.3).  Wg(3)(3.3).  Wg(5)(5.3) 

53.1 

2.6 

1.8 

knots  10  ^ft  ^ 

*s<^)24 

Wg(L)(3.3).  Wg(L)(6.2).  Wg(2)(3.5) 

35.0 

34.7 

28.0 

knots 

^Z(L)24 

Z(L)(5.3).  Z(L)(3.1).  AZ(L)(7.3) 

38.3 

6623 

5201 

ft 

-3  ,  -1 

AS.. 

b24 

Sj^(5.3).  Wg(L)(6,2),  Wg(l)(3.1) 

61.8 

2.6 

1.6 

knots  10  ft 

.  -3..-1 

AS 

a  24 

S^(5,3).Wg(2)(3.3).Sj^(7.1). 
Wg(l)(3.3).  Z(L)(3.1) 

58.6 

2.8 

1.8 

knots  10  ft 
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TABLE  VI  (Continued) 


(c)  Category  1C  (366  cases) 


Predictand 

Predictors  in  order  of  selection 

% 

red. 

Std. 

dev. 

rms 

error 

Units 

W  (L)(5,3),W  (L)(3,3),S  (5,3), 
s  s  a 

Wg(L)(4,2) 

43.2 

29.6 

22.3 

knots 

iZ(L)j2 

Z(L)(5,3),Z(L)(3,3),Z(L)(7,3), 
W^(L)(5,  1) 

35.7 

5373 

4309 

ft 

-3  -1 

Sj^(5,3),  W^(L)(4,2) 

41.6 

2.3 

1.8 

knots  10  ft 

-3  -1 

'^al2 

S  (5,3),  W  (3)(3,3) 

cl  o 

48.4 

2.5 

1.8 

knots  10  ft 

"s<'->24 

W  (1)(5,3),W  (3)(1,3),  W  (L)(6,4), 
s  s  s 

W  (3)(3,3),W  (5)(3,5),W  (2)(3,1) 

o  o  S 

37.9 

32.4 

25.5 

knots 

i2(L)24 

Z(L)(5,3),  Z(L)(7,3),  Z(L)(3,3) 

33.6 

6726 

5482 

ft 

-3  ,  -1 

^b24 

Sb(5,3),VVs(L)(6,4),W^(3)(l,3) 

55.1 

2.4 

1.6 

knots  10  ft 

-3  -1 

^a24 

S  (5,3),W  (L)(2,4),  AW  (L)(3,5) 
as  s 

52.9 

2.6 

1.8 

knots  10  ft 

(d)  Category  ID  (336  cases) 


Predictand 

Predictors  in  order  of  selection 

% 

red. 

std. 

dev. 

rms 

error 

Units 

WJ2)(5,3),W  (L)(4,2),AW  (L)(3,3), 
s  s  s 

W  (L)(6,4) 

o 

47.2 

30.6 

22.2 

knots 

AZ(L)^2 

Z(L)(5,3),Z(L)(3,3),  W^(3)(7,3) 

46.2 

7563 

5548 

ft 

-3  -1 

12 

8^(5,  3),  W^(L)(4,2),AW^(L)(3,3), 
Ws(2)(5,3) 

41.3 

2.0 

1.6 

knots  10  ft 

^"al2 

8^(5, 3),  Wg(L)(4,  2),  AWg(L)(3, 3), 
Wg(3)(5,3) 

30.2 

2.2 

1.8 

knots  lO'^ft"^ 

Ws(L)24 

Wg(2)(3, 3),  Wg(L)(5, 1),  AWg(L)(3, 3), 
Wg(2)(5,3) 

32.9 

31.7 

26.0 

knots 

AZ(L)24 

Z(L)(5,3),W  (1)(3,3),  AZ(L)(5,3) 
s 

53.5 

8228 

5611 

ft 

-3  -1 

AS„„, 
b  24 

a  24 

Sj^(5,3),  AW^(L)(7,5),  W^(L)(4,  2) 

41.5 

2.2 

1.7 

knots  10  ft 

“3  -1 

8  (5,3),  W  (2)(3,3) 
d  s 

23.7 

2.2 

2.0 

knots  10  ft 
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TABLE  VI  (Continued) 


(e)  Category  IE  (343  cases) 


Predictand 

Predictors  in  order  of  selection 

% 

red. 

Std. 

dev. 

rms 

error 

Unite 

W^(L)(5.  3),  W^(L)(4.  2),  W^(L)(7,  5). 

50.2 

28.4 

20.1 

knots 

Wg(3)(3.3) 

AZ(L)^2 

Z(L)(5.3).  Z(L)(3.3).  AZ(L)(5.3) 

38.5 

6901 

5414 

ft 

-3  -1 

^bl2 

Sj^(5,3).  Wg(L)(3,3) 

40.9 

1.9 

1.5 

knots  10  ft 

-3  -1 

AS 

a  12 

S  (5.3).  S  (7.5).  W  (3)(3.3) 
a  a  D 

38.1 

1.9 

1.5 

knots  10  ft 

*.<‘->24 

W  (L)(5.3).  W  (L)(2.4).  W  (L)(6.  2). 

D  S  S 

46.6 

30.6 

22.4 

knots 

S^(7.  5) 

a2<L)24 

Z(L)(5.3).  Z(L)(3.3).  W  (1)(7.3) 
s 

49.4 

7502 

5339 

ft 

' 

-3  -1 

AS, 

b24 

Sb(5,3).  S^(7.5) 

36.3 

2.0 

1.6 

knots  10  ft 

-3  -1 

AS 

a  24 

S  (5.3).  S  (7.5).  W  (L)(l.l) 

dL  di  B 

42.7 

2.0 

1.5 

knots  10  ft 

An  interesting  variation  is  found  for  the  most  frequent  location  on  the  predictor 
grid  of  the  selected  variables  when  category  ID  U:wo  grid  intervals  to  the  left  of  the  jet 
core)  is  compared  with  Category  IE  (two  grid  intervals  to  the  right  of  the  jet  core).  The 
second  of  the  two  numbers  in  parentheses  given  with  each  selected  variable  indicates 
whether  the  network  grid  point  is  on  the  upper  portion  of  the  grid  (numbers  4  and  5)  or 
the  lower  portion  (numbers  1  and  2).  For  Category  ID,  more  predictors  are  selected 
from  the  lower  portion  of  the  grid,  while  for  Category  IE,  more  predictors  are  selected  from 
the  upper  portion.  Thus,  for  both  categories,  the  selected  variables  are  located  on  that 
part  of  the  predictor  network  grid  that  lies  closest  to  the  jet-stream  core.  Changes  which 
take  place  in  this  region  are  normally  most  rapid  and  will  have  the  greatest  overall 
influence  on  the  future  patterns  of  the  LMW  predictands.  The  above-described  selection 
preference,  therefore,  seems  quite  logical  from  the  meteorological  viewpoint. 
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6.  Subtropical  Jet  Stream— Dependent  Data  Results 

The  dependent  sample  size  was  considerably  smaller  for  the  subtropical  jet- 
stream  categories  than  for  the  polar  categories.  Table  VII  lists  the  predictors  selected 
for  each  of  the  eight  predictands  for  each  grid-point  category,  and  the  total  percent 
reduction  of  variance  (%  red.)  of  the  predictand  for  the  predictors  selected.  Examina¬ 
tion  of  the  selected  predictors  leads  to  several  general  comments  concerning  the  selection. 
First,  as  was  the  case  with  the  polar  jet  stream,  the  first  predictor  usually  selected  for 
both  12-  and  24-hr  prognoses  was  the  initial-state  value  at  the  predictand  point  (5,3)  of 
the  LMW  parameter  being  predicted.  A  second  comment,  however,  points  to  a  significant 
difference  between  the  predictors  selected  for  the  two  jet-stream  types.  A  much  greater 
number  of  12-hr  change  predictors  were  selected  for  the  subtropical  categories.  In  many 
instances  the  change  parameter  is  the  second  predictor  selected,  indicating  the  importance 
of  the  predictor  contribution  to  the  percent  reduction  of  the  predictand.  The  fact  that  12-hr 
change  information  is  more  useful,  in  predicting  LMW  parameters  associated  with  the  sub¬ 
tropical  jet  stream,  is  not  surprising  when  one  considers  that  the  movement  of  these 
systems  is  much  more  regular  than  the  polar  jet-stream  motion  and,  therefore,  previous 
12-hr  changes  in  predictor  fields  should  have  a  higher  linear  correlation  with  future  12- 
and  24-hr  changes  of  the  predictand. 

As  would  be  expected,  aside  from  the  predictand  point  (5,  3),  the  selected  predictors 
are  most  frequently  located  at  a  point  upstream  from  the  predictand  point  for  all  categories 
except  2A.  Because  grid  points  in  this  category  fall  close  to  the  center  of  the  initial- 
state  location  of  the  wind  speed  maximum,  it  might  be  suggested  that  with  a  symmetrical 
shaped  subtropical  maximum,  predictor  information  downstream  would  be  similar  to  that 
upstream,  and,  hence,  as  likely  to  be  selected.  This  situation  would  be  much  less  likely 
to  occur  in  the  case  of  the  polar  jet  stream  where  more  rapid  changes  in  jet-stream 
intensity,  shape,  and  location  occur. 

7.  Independent-data  Testing 

The  equations  (see  the  Appendix)  obtained  with  the  selected  variables  described  in 
Sections  5  and  6,  were  tested  on  independent  data.  The  independent  data  were  obtained 
from  14  observation  times  within  the  intervals  12Z  December  16— 12Z  December  19, 

1963  and  OOZ  January  17— OOZ  January  20,  1964.  The  root-mean -square  (rms)  errors 
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TABLE  VII 

SELECTED  PREDICTORS  FOR  SUBTROPICAL  JET-STREAM  CATEGORIES 

(Results  on  dependent  data) 


(a)  Category  2A  (104  cases) 


Predictand 

Predictors  in  order  of  selection 

% 

red. 

Std. 

dev. 

rms 

error 

Units 

W  (L)(5,3).  W  (3)(5,3) 
s  s 

47.0 

32.4 

23.6 

knots 

AZ(L)^2 

Z(L)(5,  3),  W  (5)(3, 3),  W  (1)(1, 3), 

o  b 

W  (L)(7,  3),  AZ(L)(1, 3),  AW  (L)(7,  1) 

S 

68.9 

4711 

2629 

ft 

-3  -1 

'""b  12 

S^(5.3) 

50.4 

3.1 

2.2 

knots  10  ft 

-3  -1 

^"al2 

S^(5,3),W^(3)(5,3),ASj^(5,3) 

54.2 

3.4 

2.3 

knots  10  ft 

W  (L)(6,  2) 

b 

26.3 

33.4 

28.7 

knots 

AZ(L)24 

Z(L)(5,3),W  (3)(3,3),  S  (7,3) 
s  a 

51.1 

5720 

4000 

ft 

,  -3  ,  -1 

b  24 

Sj^(5.3),  AW^(L)(5,3) 

51.7 

3.4 

2.4 

knots  10  ft 

-3  -1 

a  24 

S  (5,3) 
a 

41.7 

3.7 

2.8 

knots  10  ft 

(b)  Category  2B  (90  cases) 


Predictand 

Predictors  in  order  of  selection 

%  . 

red. 

Std. 

dev. 

rms 

error 

Units 

^<^>12 

Wg(L)(5,3),Wg(L)(4,2) 

60.7 

30.9 

19.4 

knots 

AZ(L)42 

Z(L)(5,3),W  (L)(2,2) 

s 

38.8 

4458 

3487 

ft 

,  ’3  ,  -1 

^bl2 

ASJ5, 3),  AW  (L)(l,  3),  AW  (L)(7,  5), 
OS  s 

58.5 

2.3 

1.5 

knots  10  ft 

S  (5,3),  W  (L)(4,2) 
b  s 

-3  -1 

a  12 

S  (5,3),  W  (2)(5,3) 

a  b 

33.6 

1.7 

1.4 

knots  10  ft 

^s<^)24 

Wg(L)(4,  2),  Wg(L)(5,  3),  ASj^(3,  1),  S^(3, 3) 

60.0 

28.8 

18.2 

knots 

AZ(L)24 

Z(L)(5,3),  ASj^(7,5),  ASj^(3,  1), 

49.2 

4550 

3243 

ft 

S  (7,3),  S  (5,3) 

a  a 

-3  -1 

'^^b24 

Sjj(5,3),  SJ3,1) 

52.5 

2.3 

1.6 

knots  10  ft 

.1,. 

-3  -1 

^^a24 

S^(5,3),Wg(L)(4,2) 

28.8 

1.8 

1.5 

knots  10  ft 
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TABLE  VII  (Continued) 


(c)  Category  2C  (104  cases) 


Predictand 

Predictors  in  order  of  selection 

% 

red. 

Std. 

dev. 

rms 

error 

Units 

W^(LH5,3),  W^(L)(5,1) 

58.7 

31.3 

20.1 

knots 

AZ(L)^2 

AZ(LH5,3),W^(5)(3,3),  Z(L)(5,3) 

45.7 

5121 

3773 

ft 

AW^(L)(3,5) 

-3  ,  -1 

12 

S^(5,3),  AS^(7,3) 

40.6 

2.7 

2.1 

knots  10  ft 

“3  .  -1 

AS 

S  (5,3),  S  (5,1) 

43.8 

2.5 

1.9 

knots  10  ft 

a  12 

a  a 

^<^>24 

W^(2)(5,  3),  W^(L)(5, 1),  W^(L)(3,  5), 

46.3 

32.7 

24.0 

knots 

S,(5,l) 

AZ(L)24 

Z(L)(5,3) 

26.6 

5024 

4308 

ft 

-3  -1 

^^24 

S^(5,3),W^(2)(5,3),  AS^(1,3), 

55.9 

2.8 

1.8 

knots  10  ft 

AZ(L)(3,3),AS  (1,1),W  (3)(3,1) 
a  s 

-3  -1 

^^a24 

S  (5,3) 
a 

42.2 

2.6 

2.0 

knots  10  ft 

(d)  Category  2D  (138  cases) 


Predictand 

Predictors  in  order  of  selection 

% 

red. 

Std. 

dev. 

rms 

error 

Units 

^<^>12 

W  (L)(5,  3),  W  (2)(3,  3),  W  (5)(5, 3) 

o  S  S 

59.3 

30.3 

19.3 

knots 

AZ(L)^2 

Z(L)(5,3),  Z(L)(3,3) 

20.8 

4788 

4261 

ft 

-3  -1 

^"bl2 

S^(5,3),  S^(3,3) 

31.7 

1.9 

1.6 

knots  10  ft 

-3  -1 

^"al2 

S  (5,3),  W  (L)(3,3) 
a  b 

46.1 

2.0 

1.4 

knots  10  ft 

'^s(^)24 

W  (L)(4,2),  W  (2)(3,3) 

O  b 

38.4 

34.2 

26.8 

knots 

AZ(L)24 

Z(L)(5,3),  W  (5)(1,3) 
s 

33.0 

5770 

4757 

ft 

~3  -1 

'^®b24 

S^^(5,3),  W^(L)(4,2) 

37.3 

2.3 

1.9 

knots  10  ft 

"3  "1 

^"a24 

S  (5,3),  W  (L)(4,2) 
a  s 

46.9 

2.4 

1.7 

knots  10  ft 
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TABLE  VII  (Continued) 


(e)  Category  2E  (98  cases) 


Predictand 

Predictors  in  order  of  selection 

% 

red. 

Std. 

dev. 

rms 

error 

Units 

W  (3)(5,  3),  W  (1)(7, 3),  w  (5)(3, 3) 

S  b  D 

50.6 

26.3 

18.5 

knots 

AZ(L)^2 

Z(L)(5,3),W^(3)(5,3),Z(L)(3,3), 

s 

W  (5)(5,3),  W  (L)(3,3) 

S  D 

57.6 

5046 

3284 

ft 

,  -3  -1 

^bl2 

Sjj(5,3) 

29.2 

2.4 

2.0 

knots  10  ft 

-3.-1 

^"al2 

S  (5,3).  W  (5)(7.3) 

d.  b 

57.8 

3.1 

2.0 

knots  10  ft 

^<^>24 

W  (3)(5.3).  W  (L)(6.4),  AZ(L)(3.3) 
s  s 

53.1 

26.1 

17.9 

knots 

AZ(L)24 

Z(L)(5.3) 

46.8 

5524 

4027 

ft 

-3  -1 

“^^5  24 

S  (5.  3).  AS  (3, 3).  AW  (L)(3,3) 

D  a  s 

52.7 

2.2 

1.5 

knots  10  ft 

.  -3  -1 

24 

S  (5.3).  AS  (5.3) 
a  a 

51.5 

2.0 

1.4 

knots  10  ft 

of  the  equations  for  12-  and  24-hour  prognoses  are  compared  with  persistence  (forecast 

* 

of  no  change)  for  all  categories  of  the  polar  jet  stream  in  Table  VIII.  The  independent 
sample  size  for  the  five  grid-point  categorie's  varies  between  82  and  85  cases.  The  grid 
points  in  the  independent  data  at  which  the  prognoses  are  made  are  selected  for  each 
category  with  the  same  objective  procedure  that  was  used  for  the  dependent  data  sample. 

The  prognoses  obtained  from  the  derived  equations  are  superior  to  persistence 
(lower  rms  error)  for  all  categories,  both  time  intervals,  and  each  predictand.  In  most 
instances,  the  rms  error  is  significantly  lower  than  persistence.  The  most  outstanding 
example  is  found  in  the  Category  lA,  24-hr  prognoses  of  wind  speed,  shear  above,  and 
shear  below;  here,  the  rms  errors  of  the  equations  are  all  less  than  half  that  of  persistence. 
The  persistence  rms  errors  are  highest  for  this  category  and  forecast  interval.  This  is  to 
be  expected,  as  the  position  and  intensity  of  the  wind  speed  maximum  can  change  drastically 
in  a  24-hr  period. 

If  the  independent-data  equation  rms  errors  shown  in  Table  VIII  are  compared  to 
the  dependent-data  rms  errors  it  is  seen  that  the  two  sets  of  error  statistics  are  very 
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TABLE  Vm 

Rms  ERRORS  ON  INDEPENDENT  DATA  FOR 
POLAR  JET-STREAM  CATEGORIES 


(a)  Category  lA  (84  cases) 


Predictand 

Prediction 

rms  error 

Persistence 

rms  error 

Units 

26.2 

48.4 

knots 

AZ(L)^2 

4785 

5624 

ft 

,  -3  -1 

^^bl2 

1.7 

3.1 

knots  10  ft 

-3  -1 

^\l2 

1.9 

3.6 

knots  10  ft 

*.<^>24 

30.4 

61.9 

knots 

AZ(L)24 

5101 

7041 

ft 

-3  -1 

b  24 

1.4 

3.6 

knots  10  ft 

-3  -1 

a  24 

2.0 

4.1 

knots  10  ft 

(b)  Category  IB  (84  cases) 


Predictand 

Prediction 

rms  error 

Persistence 

rms  error 

Units 

24.5 

34.2 

knots 

AZ(L)^2 

3731 

4974 

ft 

-3  ,  -1 

^^bl2 

1.6 

2.5 

knots  10  ft 

-3  .  -1 

^al2 

1.7 

2.8 

knots  10  ft 

26.9 

42.9 

knots 

AZ(L)24 

4928 

6064 

ft 

-3  -1 

'^b24 

1.6 

2.7 

knots  10  ft 

-3  -1 

^a24 

1.9 

3.3 

knots  10  ft 
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TABLE  Vni  (Continued) 


(c)  Category  1C  (85  cases) 


Predictand 

Prediction 

rms  error 

Persistence 

rms  error 

Units 

^s(^)l2 

23.9 

28.0 

knots 

AZ(L)^2 

4770 

5520 

ft 

-3.  -1 

^^12 

1.7 

2.2 

knots  10  ft 

-3  ^ 

AS 

a  12 

2.1 

2.7 

knots  10  ft 

^(^>24 

32.3 

38.7 

knots 

az(L)2^ 

4218 

5380 

ft 

--3  -1 

AS, 
b  24 

2.2 

2.6 

knots  10  ft 

-3  -1 

a  24 

2.2 

3.2 

knots  10  ft 

(d)  Category  ID  (83  cases) 


Predictand 

Prediction 

rms  error 

Persistence 

rms  error 

Units 

^s(^)l2 

24.0 

30.1 

knots 

AZ(L)^2 

5727 

7837 

ft 

-3  -1 

^‘bl2 

1.5 

2.2 

knots  10  ft 

-3  -1 

^al2 

1.4 

1.8 

knots  10  ft 

33.6 

40.0 

knots 

AZ(L)24 

5714 

8022 

ft 

-3  -1 

^^24 

2.2 

2.7 

knots  10  ft 

-3  -1 

^a24 

2.2 

2.4 

knots  10  ft 
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TABLE  VIII  (Continued) 


(e)  Category  IE  (82  cases) 


Predictand 

Prediction 

rms  error 

Persistence 

rms  error 

Units 

17,6 

19.0 

knots 

AZ(L)^2 

5353 

5604 

ft 

-3  -1 

1.3 

1.7 

knots  10  ft 

,  -3  ,  -1 

^^al2 

1.5 

1,7 

knots  10  ft 

23.4 

25.3 

knots 

AZ(L)24 

5059 

5323 

ft 

“3  ,  -1 

AS, 
b  24 

1.6 

2.0 

knots  10  ft 

-3  -1 

AS 

a  24 

1.5 

1.9 

knots  10  ft 
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similar,  indicating  the  stability  of  the  relationships  on  the  independent-data  sample 
tested. 

The  rms  errors  on  independent  data,  obtained  with  the  12-  and  24-hr  subtropical 
jet-stream  prediction  equations  for  all  categories,  are  given  in  Table  IX.  The  sample 
size  varies  between  17  and  28  cases  and,  because  of  the  limited  size,  the  results  must  be 
interpreted  cautiously.  However,  the  results  appear  consistent  with  a  priori  meteorological 
reasoning  and  will  be  presented  in  this  light. 

The  persistence  rms  errors  for  24-hr  prediction  of  wind  speed,  shear  above,  and 
shear  below  are  again  highest  in  the  wind  speed  maxima  category  (2A).  The  24-hr  per¬ 
sistence  rms  errors  for  wind  speed  and  height  are,  for  most  categories,  lower  for  the 
subtropical  jet  stream  than  those  obtained  for  the  polar  jet  stream.  This  result  is  con¬ 
sistent  with  the  more  regular  behavior  and  persistence  of  the  subtropical  jet  stream.  How¬ 
ever,  no  consistent  difference  can  be  noted  between  the  rms  errors  for  the  two  types  of 
jet  streams,  when  considering  shear  above  and  below  persistence  forecasts. 

The  rms  errors  obtained  from  the  prediction  equations  are  lower  than  persistence 
with  only  one  exception.  The  rms  error  for  the  12-hr  forecast  of  shear  above  for 
Category  2B  is  the  same  as  persistence.  A  comparison  of  the  independent-data  rms 
errors  given  in  Table  IX  with  the  dependent-data  rms  errors  in  Table  VII  shows  that,  in 
general,  the  rms  errors  obtained  with  independent  data  are  similar,  although  much  higher 
or  lower  errors  do  occur  in  some  categories.  These  fluctuations  are  very  likely  a  result 
of  the  small  sample  size  of  the  independent  data. 

A  final  series  of  experiments  was  conducted  with  the  independent  data  to  consider 
the  usefulness  of  the  categorization  that  was  employed  in  this  feasibility  study.  Specifically, 
the  three  features  investigated  were:  (1)  separation  into  polar  and  subtropical  jet-stream 
classes,  (2)  the  use  of  exclusive  categories  to  the  right  and  the  left  of  the  jet  core  and 
(3)  the  breakdown  of  categories  along  the  core. 

Table  X  gives  the  rms  errors  that  occur  when  the  equation  developed  for  one 
category  is  tested  on  the  independent  data  cases  of  a  different  category.  A  comparison 
with  the  rms  error  statistics  shown  in  Tables  VIII  and  IX  is  required  in  the  interpretation 
of  the  results. 

Parts  (a)  and  (b)  of  Table  X  provide  information  concerning  items  (1)  and  (2)  above. 
The  equations  developed  for  Category  ID  (two  grid  intervals  to  the  left  of  the  polar  jet- 
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TABLE  EX 

Rms  ERRORS  ON  INDEPENDENT  DATA  FOR 
SUBTROPICAL  JET-STREAM  CATEGORIES 


(a)  Category  2A  (20  cases) 


Predictand 

Prediction 

rms  error 

Persistence 

rms  error 

Units 

24.9 

41.0 

knots 

AZ(L)^2 

3187 

4897 

ft 

-3 ,  -1 

^"bl2 

1.7 

3.2 

knots  10  ft 

-3.  -1 

^al2 

1.9 

3.3 

knots  10  ft 

^(^>24 

27.9 

40.3 

knots 

AZ(L)24 

3353 

4856 

ft 

,  -3  ,  -1 

AS 

b24 

2.3 

3.4 

knots  10  ft 

-3  ,  -1 

a  24 

2.6 

4.1 

knots  10  ft 

(b)  Category  2B  (20  cases) 


Predictand 

Prediction 

rms  error 

Persistence 

rms  error 

Units 

28.6 

32.8 

knots 

AZ(L)^2 

3100 

3668 

ft 

-3 ,  -1 

^bl2 

2.7 

3.2 

knots  10  ft 

“3  ,  -1 

^^12 

2.4 

2.4 

knots  10  ft 

^("^>24 

21.8 

29.3 

knots 

AZ(L)24 

3327 

4853 

ft 

.  -3.-1 

^^24 

2.1 

2.9 

knots  10  ft 

-3 .  -1 

^"a24 

1.5 

1.9 

knots  10  ft 
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TABLE  EX  (Continued) 


(c)  Category  2C  (18  cases) 


Predictand 

Prediction 

rms  error 

Persistence 

rms  error 

Units 

28.2 

33.4 

knots 

AZ(L)^2 

4180 

5780 

ft 

-3 ,  -1 

^Sbl2 

1.8 

2.2 

knots  10  ft 

-3  ,  -1 

^Sal2 

1.9 

2.5 

knots  10  ft 

27.5 

37.6 

knots 

AZ(L)24 

4770 

6800 

ft 

-3 ,  -1 

AS 

b24 

2.0 

3.3 

knots  10  ft 

-3 ,  -1 

AS 

a  24 

1.9 

2.8 

knots  10  ft 

(d)  Category  2D  (28  cases) 


Predictand 

Prediction 

rms  error 

Persistence 

rms  error 

Units 

20.3 

25.1 

knots 

AZ(L)^2 

4616 

5904 

ft 

-3  -1 

^bi2 

1.4 

2.0 

knots  10  ft 

-3  -1 

AS 

a  12 

1.5 

2.3 

knots  10  ft 

^(^>24 

25.1 

31.4 

knots 

AZ(L)24 

4570 

6508 

ft 

-3  -1 

^b24 

1.9 

2.5 

knots  10  ft 

-3  ,  -1 

^a24 

1.6 

2.5 

knots  10  ft 
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TABLE  IX  (Continued) 


(e)  Category  2E  (17  cases) 


Predictand 

Prediction 

rms  error 

Persistence 

rms  error 

Units 

21.5 

22.3 

knots 

AZ(L)^2 

5022 

6933 

ft 

-3  .  -1 

^^bl2 

1.6 

2.2 

knots  10  ft 

-3  ,  -1 

AS 

a  12 

1.5 

2.6 

knots  10  ft 

23.9 

28.1 

knots 

AZ(L)24 

2992 

4011 

ft 

-3 

^b24 

1.7 

3.0 

knots  10  ft 

-3 

'^a24 

1.9 

3.0 

knots  10  ft 
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TABLE  X 

RMS  ERRORS  ON  INDEPENDENT  DATA  FOR 
COMPARISON  OF  CATEGORY  EQUATIONS 

(a)  Category  lA  Data  (84  cases) 


Rms  error 

Predictand 

Category  ID 
equation 

Category  2A 
equation 

Units 

30.5 

32.5 

knots 

AZ(L)^2 

5121 

5703 

ft 

-3  .  -1 

^"bl2 

1.9 

1.8 

knots  10  ft 

,  -3  ,  -1 

^^12 

2.1 

2.4 

knots  10  ft 

34.9 

38.8 

knots 

AZ(L)24 

5235 

5412 

ft 

-3  -1 

^"b24 

1.4 

1.6 

knots  10  ft 

“3  -1 

24 

2.2 

2.1 

knots  10  ft 

(b)  Category  ID  Data  (83  cases) 


Rms  error 

Predictand 

Category  LA 
equation 

Category  IE 
equation 

Units 

35.4 

27.4 

knots 

AZ(L)^2 

6293 

6060 

ft 

^Sbl2 

1.6 

1.6 

•“3  -1 

knots  10  ft 

^Sal2 

1.9 

1.6 

-3  _i 

knots  10  ft 

^<"^^24 

39.8 

35.0 

knots 

AZ(L)24 

6144 

6501 

ft 

b  24 

2.3 

2.3 

-3  -1 

knots  10  ft 

a  24 

2.4 

2.2 

knots  10  ^  ft  ^ 
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TABLE  X  (Continued) 


(c)  Category  IB  Data  (84  cases)  and 
Category  1C  Data  (85  cases)* 


Rms  error 

Predictand 

Category  1C 
equation 

Category  IB 
equation 

Units 

22.6 

25.4 

knots 

AZ(L)i2 

3537 

5106 

ft 

-3  -1 

^^bl2 

1.5 

1.7 

knots  10  ft 

-3  -1 

^"al2 

1.7 

2.0 

knots  10  ft 

25.7 

33.3 

knots 

AZ(L)24 

5196 

4311 

ft 

-3  -1 

"^^b24 

1.6 

2.2 

knots  10  ft 

-3  -l' 

^\24 

1.9 

2.4 

knots  10  ft 

(d)  Combined  Data  of  Categories  2A,  2B,  and  2C  (56  cases) 


Rms  error 

Predictand 

Category  2A 
equation 

Category  2B 
equation 

Category  2C 
equation 

Units 

29.1 

29.6 

28.0 

knots 

AZ(L)i2 

3954 

3628 

3882 

ft 

-3  -1 

^^bl2 

2.1 

2.4 

2.4 

knots  10  ft 

-3  -1 

^"al2 

2.0 

2.5 

2.2 

knots  10  ft 

'^s(^)24 

26.6 

26.3 

28.4 

knots 

AZ(L)24 

3777 

4104 

4178 

ft 

-3  -1 

^®b24 

2.2 

2.2 

2.5 

knots  10  ft 

-3  -1 

^"a24 

2.0 

2.2 

2.0 

knots  10  ft 

*Category  IB  equation  tested  on  Category  1C  data,  and  vice  versa. 
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stream  core)  and  Category  2A  (subtropical  wind  speed  maximum)  are  tested  on  inde¬ 
pendent  data  for  Category  lA  (polar  wind  speed  maximum).  The  rms  errors  for  these 
two  equations  [Table  X(a)]  are  consistently  higher  for  all  predictands  than  the  rms  errors 
for  the  Category  lA  equation  [Table  Vlll(a)].  For  example,  the  rms  errors  for  24-hr  prog¬ 
noses  of  wind  speed  are  30.4  for  the  Category  lA  equation,  34,9  for  the  Category  ID 
equation,  and  38.8  for  the  Category  2A  equation.  Thus,  the  equations  developed  for  the 
wind-speed  maximum  category  of  the  subtropical  jet  and  for  a  category  off  the  polar 
jet  core  yielded  inferior  results  when  tested  on  independent  data  for  the  polar  jet- 
stream  wind- speed  maximum  category. 

The  rms  errors  obtained  when  the  Category  lA  equation  and  Category  IE  equation 
(developed  for  two  grid  intervals  to  the  right  of  the  polar  jet-stream  core)  are  tested  on 
Category  ID  independent  data  are  given  in  Table  X(b),  A  comparison  with  the  rms  errors 
obtained  with  the  Category  ID  equation  [Table  Vin(d)]  shows  that  higher  rms  errors  resulted 
for  all  predictands  using  both  the  Category  lA  and  Category  IE  equations.  Taking  the  24-hr 
wind  speed  as  an  example  again,  the  rms  errors  are  33.6  for  the  Category  ID  equation, 

39.8  for  the  Category  LA  equation,  and  35,0  for  the  Category  IE  equation.  It  is  seen, 
therefore,  that  the  equations  (Category  lA  equation)  developed  for  a  category  on  the  core 
yield  increased  rms  errors  when  applied  to  an  area  removed  from  the  core.  Furthermore, 
when  the  equations  developed  for  the  category  to  the  right  of  the  core  are  applied  to  data 
to  the  left  of  the  core,  higher  rms  errors  result.  The  above  results  indicate  that  the 
stratification  of  the  data  into  polar  and  subtropical  systems,  each  being  further  separated 
into  grid  points  on  the  core  or  to  the  right  or  left  of  the  core,  was  worthwhile, 

Rms  errors  that  occur  when  the  Category  1C  equation  (two  grid  points  downstream 
from  the  wind  speed  maximum)  is  applied  to  Category  IB  independent  data,  and  vice 
versa,  are  given  in  Table  X(c),  If  the  classification  procedure  used  along  the  core  has 
merit,  considerably  higher  rms  errors  should  result  when  the  equation  developed  for  the 
category  upstream  from  the  wind  speed  maximum  is  applied  to  independent  data  down¬ 
stream,  A  comparison  with  Table  Vlll(b)  and  (c)  shows  that  this  is  not  the  case.  The  rms 
errors  differ  very  little,  and  for  some  predictands  are  actually  lower.  This  result  may 
be  due  to  the  fact  that  wind  speed  and  shear  variations  along  the  core  are  much  smaller 
than  the  variations  which  occur  to  the  right  and  left  of  the  core.  The  obvious  conclusion 
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to  be  drawn  from  the  above  is  that  a  single  category  of  grid  points  located  on  the  jet- 
stream  core  should  replace  the  three  categories  used.  The  results  shown  in  Table  X{d) 
further  substantiate  this  conclusion.  The  sets  of  equations  for  the  three  subtropical 
core  categories  are  applied  to  all  grid  points  along  the  subtropical  jet  core,  with  little 
difference  noted  in  the  rms  errors. 

An  example  of  the  predictions  obtained  from  the  LMW  Predictor  Network  Technique 

(PNT)  is  given  in  Fig.  4(a)— (h).  The  initial-state  distribution  of  W  (L)  on  OOZ  December 

s 

17,  1963  is  given  in  Fig.  4(a).  The  subsequent  12-  and  24-hr  distributions  (analyses 
obtained  from  observations)  are  shown  in  Fig.  4(b)  and  (d).  The  corresponding  12-  and 
24-hr  prognoses  are  given  in  Fig.  4(c)  and  (e).  The  isopleths  for  these  maps  and 
Fig.  4(f),  (g)  and  (h)  (24-hr  prognoses  of  height,  shear  below,  and  shear  above)  are  all  drawn 
subjectively.  The  isopleths  were,  however,  based  on  grid-point  values  of  the  parameters 
that  were  obtained  objectively  from  the  prediction  equations  or,  in  the  case  of  the  LMW 
wind-speed  analyses,  from  a  SAT  analysis  of  the  observations.  The  jet-stream  cores 
shown  on  the  LMW  wind  speed  analyses  and  prognoses  are  drawn  in  subjectively,  based 
on  the  grid-point  values  of  wind  speed  and  the  isotach  analysis.  Both  the  predicted  and 
analyzed  jet  cores  are  shown  on  the  maps  depicting  the  24-hr  predictions  of  LMW  height, 
shear  above  and  shear  below. 

An  examination  of  the  initial-state  LMW  wind  speed  map  shows  that  three  major, 
intense  jet-stream  systems  with  wind  maxima  greater  than  125  knots  are  to  be  found. 

The  system  over  the  southeastern  North  Pacific  Ocean  curves  southeastward  toward  the 
equator.  A  second  strong  system  defines  an  upper-level  ridge  centered  just  off  the 
Pacific  coast  of  the  U.S.  and  southwest  Canada.  The  isotach  maximum  is  moving  across 
the  top  of  the  ridge  and  is  located  at  the  Canadian  border.  The  third  system  is  a  major 
east-west  jet  stream  extending  from  the  eastern  United  States  across  the  Atlantic  just 
south  of  40°N  latitude.  An  intense  wind-speed  maximum  of  200  knots  is  noted  at  50°W 
longitude. 

The  major  changes  that  occur  after  24  hours  are:  (1)  the  jet  stream  over  the 
southeastern  North  Pacific  Ocean  intensifies  upstream  at  the  edge  of  the  map  and  begins 
to  curve  slightly  north  of  east,  tending  to  be  linked  with  the  system  over  the  west  coast 
of  the  United  States,  (2)  the  jet  stream  near  the  west  coast  of  the  U.S.  is  displaced  north- 
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ward  at  the  crest  of  the  ridge,  but  the  jet  maximum  plunges  into  the  northern  plains, 
and  (3)  the  jet  stream  in  the  Atlantic  maintains  position  but  the  intensity  is  somewhat 
diminished. 

The  corresponding  24-hr  prognoses  (1)  correctly  forecast  the  change  in  curvature  of 
the  Pacific  Ocean  jet  stream  and  the  tendency  to  ‘link”  with  the  west  coast  jet-stream  system, 
but  do  not  indicate  the  increased  intensity  upstream,  (2)  correctly  indicate  the  departure  of 
the  Pacific  Coast  wind- speed  maximum  south-eastward,  but  do  not  displace  the  jet-stream 
northward  and  (3)  correctly  maintains  the  position  and  reduces  the  intensity  of  the  Atlantic 
jet  stream  system,  but  does  not  indicate  the  jet  maximum  in  the  extreme  eastern  Atlantic 
Ocean. 

In  general  then,  it  can  be  said  that  the  prediction  equations  for  the  LMW  wind  speeds 
yield  quite  reasonable  prognoses.  A  weakness  in  the  prediction  technique  is,  however, 
clearly  illustrated  by  these  series  of  maps.  A  need  exists  to  make  predictions  further 
downstream  from  the  initial-state  location  of  the  jet-stream  wind-speed  maximum  than 
was  undertaken  with  the  technique  developed  in  this  feasibility  study  (no  predictions  were 
made  over  most  of  the  western  and  central  United  States  because  no  predictand  points 
were  located  in  those  areas;  this  resulted  in  the  LMW  parameters  being  undefined  in  this 
region  of  the  prognoses). 

Figure  4(f),  (g)  and  (h)  shows  the  24-hr  predictions  of  the  LMW  shear  below,  shear 
above  and  height  of  the  LMW  with  the  actual  and  forecast  jet  cores  superimposed  on  the 
analyses.  Generally,  the  most  intense  shears  above  the  LMW  are  located  close  to  or 
just  south  of  the  jet  core,  while  the  most  intense  shears  below  the  LMW  are  located  close 
to  or  just  north  of  the  core.  The  LMW  height  prognoses  are  shown  only  within  the  area 
defined  by  the  predicted  75- knot  isotach. 

To  summarize.  Fig.  4  shows  that  reasonable  forecasts  are  obtained  for  a  particular 
series  of  maps.  The  series  of  maps  also  helps  to  emphasize  the  point  that  additional 
categories  of  predictand  points  over  areas  where  there  are  no  jets  at  initial  time  (t^) 
would  be  required  before  developing  an  operationally-useful  set  of  prediction  equations. 
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Fig.  4(a).  W  (L)  analysis,  OOZ  December  17,  1963. 
s 
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Fig.  4(d).  W  (L)  analysis,  OOZ  December  18,  1963. 
s 


Fig.  4(e).  W  (L)  24-hr  prognosis  valid  OOZ  December  18,  1963. 
s 
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Fig.  4(f). 


S,  24-hr  prognosis  valid  OOZ  December  18,  1963. 
b 


Fig.  4(g).  S  24-hr  prognosis  valid  OOZ  December  18,  1963. 

3. 
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Fig.  4(h).  Z(L)  24-hr  prognosis  valid  OOZ  December  18,  1963. 
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SECTION  V 


RESULTS  USING  THE  MODELING  PREDICTION  TECHNIQUE 

The  LMW  Predictor  Network  Technique  was  developed  from  a  two-month  winter 
data  sample,  a  portion  of  which  was  set  aside  as  independent  data  for  testing  the 
technique.  Prognostic  constant -pressure-surface  charts  were  not  included  as  part 
of  this  data  sample.  Thus,  it  was  necessary  to  acquire  a  special  collection  of  numeri¬ 
cal  prognoses  for  the  LMW  Modeling  Prediction  Technique,  which  requires  using 
constant-pressure-surface  prognostic  charts  in  conjunction  with  the  diagnostic  LMW 
equations.  Third  Weather  Wing  personnel  (Offutt  Air  Force  Base,  Omaha,  Nebraska) 
supplied  us  with  a  five-day  sample  of  prognostic  charts  and  analyses  (9—14  December 
1964)  for  testing. 

Ten  sets  of  12-,  24-,  36-,  and  48-hr  constant-pressure-surface  prognoses  of 
height,  temperature,  and  wind  were  used  with  the  diagnostic  (as  opposed  to  the  prognostic 
equations  discussed  in  the  previous  section)  LMW  equations,  beginning  with  the  set 
that  was  generated  from  initial  data  at  OOZ  December  9,  1964.  Analyses  of  height, 
temperature,  and  wind  from  OOZ  December  9  to  OOZ  December  14  were  used  with 
the  diagnostic  LMW  equations  to  obtain  LMW  analyses  for  verification  purposes  (the 
analyses  were  not  “corrected”  by  station  data  and  the  successive  approximation 
analysis  technique  as  described  in  PART  I  of  this  report,  because  station  data  were 
not  on  hand).  Because  analyses  were  not  available  after  OOZ  December  14,  ten  12-hr, 
nine  24-hr,  eight  36-hr,  and  seven  48-hr  LMW  prognoses  were  verified  with  the 
analyses. 

Root-mean-square  errors  were  computed  (at  grid  points)  for  each  forecast 
LMW  parameter  for  each  observation  time,  and  overall  rms  errors  were  computed 
for  each  LMW  parameter  for  the  12-,  24-,  36-,  and  48-hr  prognosis  groups.  As  was 
expected,  the  errors  increased  with  increasing  length  of  the  forecast  interval.  Table 
XI  presents  the  overall  error  statistics. 

The  measure  of  error  by  rms  errors  for  the  prognoses  of  the  LMW  parameters 
is  not  a  particularly  good  indicator  of  their  accuracy,  because  of  the  nature  of  the 
horizontal  variations  of  the  fields  of  LMW  parameters. 
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TABLE  XI 


OVERALL  rms  ERRORS  FOR  LM\V  PROGNOSES 


Predictand 

rms  errors 

Units 

12-hr 

24-hr 

36-hr 

48-hr 

W  (L) 
s 

12.8 

15.7 

18.6 

21.0 

knots 

Z  (L) 

4598 

5289 

5924 

6491 

ft 

-3  -1 

0.98 

1.04 

1.09 

1.13 

knots  10  ft 

S 

-3  -1 

0.92 

1.01 

1.11 

1.20 

knots  10  ft 

a 

To  obtain  a  clearer  picture  of  the  usefulness,  the  prognoses  were  examined  with 
regard  to:  (1)  location  of  the  jet-stream  axes,  (2)  the  intensity  of  the  maxima  along 
the  axes,  (3)  the  extent  of  the  areas  of  wind  speeds  in  excess  of  100  knots,  and 
(4)  the  height  of  the  LMW  and  the  vertical  shears  below  and  above  the  LMV^  over  those 
areas  where  wind  speeds  were  observed  in  excess  of  100  knots.  The  evaluation  based 
on  this  type  of  examination  was  necessarily  subjective  in  nature.  Accuracy  of  the 
predictions  of  each  of  the  LMW  parameters  is  discussed  separately. 

8.  Wind  speed  at  the  LMW  [W  (L)1 

s 

The  12-hr  prognoses  of  W^(L)  were  considered  very  good,  based  on  the  first 
three  criteria  listed  above,  i.e.  (1)  the  location  of  the  jet-stream  axes  was  generally 
forecast  quite  well  with  only  occasional  maximum  errors  in  position  of  the  jet-cores 

of  100—150  miles,  (2)  errors  for  the  intensity  of  the  maxima  along  the  isotachs 

,  r  .  ,nr  ,Fcst  minus  observed . . . 

ranged  from  near  zero  to  25%  (%  error  =  p - ; - ; - 1  X  100)  with  the 

observed 

average  errors  about  10%,  and  (3)  the  extent  of  the  areas  of  wind  speeds  in  excess 
of  100  knots  was  specified  very  well  over  most  areas  of  the  hemisphere  for  the  12-hr 
prognoses. 

The  24-hr  prognoses  of  W^(L)  were  judged,  in  general,  to  be  fair-to-good  or 
good.  The  locations  of  the  jet-stream  axes  were  forecast  well  over  most  areas  of 
the  hemisphere  during  most  of  the  five-day  period,  but  over  central  and  eastern  Asia 
and  the  western  Pacific,  where  there  were  large-amplitude  trough— ridge  systems  in 
the  planetary  waves,  causing  the  jet-streams  axes  to  have  a  markedly  meridional 
character,  the  extent  of  the  meridional  air  flow  was  usually  under-forecast,  although 
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the  relative  longitudinal  positions  of  the  troughs  and  ridges  were  forecast  accurately. 
Errors  for  the  intensity  of  the  maxima  were  somewhat  higher  than  for  the  12-hr 
forecasts,  ranging  up  to  40%  occasionally,  but  averaging  15—20%. 

The  area  of  winds  in  excess  of  100  knots  was  specified  fairly  accurately  over 
most  areas  of  the  hemisphere,  with  no  apparent  systematic  error.  At  times,  there 
were  areas  where  winds  were  forecast  to  be  in  excess  of  100  knots,  but  were  not 
observed  to  be  and  the  converse  was  true  at  other  times.  In  other  words,  there  was 
no  noticeable  bias  in  the  forecasts. 

The  36-hr  and  48-hr  prognoses  of  W  (L)  showed  considerable  deterioration  from 

s 

the  12-  and  24-hr  prognoses,  although  two  prognoses  in  each  group  were  considered 

to  be  good,  with  the  remainder  only  fair  to  fair-to-good.  The  largest  errors  in  the 

forecasts  occurred  over  central  and  eastern  Asia  where  large  amplitude  planetary 

waves  were  evident,  and  over  an  area  near  Hawaii  in  the  vicinity  of  a  cut-off  closed 

cyclonic  circulation  that  was  quasi-stationary  for  half  the  five-day  period  and  then 

began  to  move  northeastward.  Over  the  North  Atlantic  Ocean,  where  a  wide  area 

of  strong  zonal  (westerly)  winds  (^100  knots)  was  observed  for  most  of  the  five-day 

period,  the  36--  and  48-hr  prognoses  were  judged  to  be  good. 

Figure  5(a)  through  (e)  is  the  observed  and  12-,  24-,  36-,  and  48-hr  W  (L) 

s 

prognoses  valid  at  OOZ  December  13.  These  prognoses  were  judged  to  be  fair-to- 
good.  Figure  6(a)  through  (e)  is  the  12-,  24-,  36-,  and  48-hr  prognoses  valid  at 
OOZ  December  14— a  series  of  prognoses  that  were  considered  to  be  generally  good- 
to -excellent. 

9.  Height  of  the  LMW  [Z(L)] 

The  areas  where  the  wind  speeds  were  observed  in  excess  of  100  knots  were 

examined  on  the  prognosis  charts  of  Z(L)  and  compared  with  the  observed  chart 

of  Z(L).  It  was  evident  from  the  comparison  that  the  errors  in  the  prognoses  were 

relatively  low — generally  less  than  2000  ft  and  rarely  more  than  5000  ft  for  the 

12-  and  24-hr  prognoses  and  somewhat  higher  (2000—4000  ft  and  only  occasionally  as 

'high  as  8000—10000  ft)  for  the  36-  and  48-hr  prognoses.  In  general,  the  prognoses 

of  Z(L)  for  all  forecast  intervals  were  considered  to  be  good.  Figure  7(a)  through 

(d)  shows  the  error  fields  for  Z(L)  prognoses  [over  areas  where  W  (L)  was  observed 

s 

as  2:100  knots]  valid  at  OOZ  December  13,  1964,  and  Fig.  8(a)  through  (d)  shows  the 
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Z(L)  error  fields  for  prognoses  valid  24  hr  later.  (These  sets  of  error  fields  correspond 

to  the  verification  times  for  the  W  (L)  prognoses  discussed  in  the  previous  section.) 

s 

It  is  evident  that  the  general  conclusion  stated  for  all  the  Z(L)  prognoses  in  the  five-day 
series  is  valid  for  these  specific  observation  times. 

10.  Shear  below  the  LMW  (S,  ) 

b 

Within  the  100-knot  wind-speed  isotach,  errors  generally  averaged  about 
1.5  knots  per  thousand  feet  for  the  12-hr  prognoses  and  close  to  2  knots  per  thousand 
feet  for  the  24-hr  prognoses  except  for  very  small  areas  where  the  errors  were 
sometimes  as  high  as  3  to  6  knots  per  thousand  feet.  It  should  be  emphasized  that  the 
areas  where  these  larger  errors  occurred  comprised  only  about  5%  of  the  areas  enclosed 
by  the  100-knot  isotachs.  For  the  36-  and  48-hr  prognoses,  the  average  errors  over 
the  areas  within  the  100-knot  isotach  increased  to  2  to  3  knots  per  thousand  feet,  but 
maximum  errors  of  4  to  6  knots  occurred  only  over  a  small  portion  (3—5%)  of  the 
total  area  outlined  by  the  observed  100-knot  isotach.  In  contrast,  the  total  rms  errors 
for  the  hemisphere  for  all  forecast  intervals  was  near  1  knot  per  thousand  feet, 
because  over  the  major  portion  of  the  hemisphere  winds  are  less  than  100  knots  at 
the  LMW,vertieal  shears  are  less,  and  it  follows  that  the  errors  are  quite  small. 

The  areas  of  maximum  errors  of  appear  to  be  highly  correlated  to  the  maximum 
errors  of  wind  speed  at  the  LMW„ 

Figure  9('a)  through  (d)  and  Fig  10.  (a)  through  (d)  present  the  error  distribution 

of  forecast  fields  of  within  the  observed  100-knot  isotach  for  the  same  time  periods 

as  for  W  (L)  and  Z(L). 
s 

11.  Shear  Above  the  LMW  (S  ) 

-  a 

The  statements  in  the  discussion  of  the  shear  below  the  LMW  apply  for  S  , 

3. 

except  that  there  were  a  few  cases  at  isolated  grid  points  where  maximum  errors 

were  8  and  9  knots  per  thousand  feet.  Figure  11(a)  through  (d)  and  Fig  12  (a)  through 

(d)  are  the  error  fields  within  the  100-knot  wind-speed  isotach  of  S  prognoses  valid 

a 

at  OOZ  December  13  and  OOZ  December  14. 
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Ill 


Fig.  5(c).  W  (L)  24-hr  prognosis  valid  OOZ  December  13,  1964. 
s 
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Fig.  5(d).  W  (L)  36-hr  prognosis  valid  OOZ  December  13,  1964. 
s 
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Fig.  6(b).  W  (L)  12-hr  prognosis  valid  OOZ  December  14,  1964. 
s 


Isotach,  knots 
Error  field,  knots 
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Fig.  6(c).  W^(L)  24-hr  prognosis  valid  OOZ  December  14,  1964. 


Fig.  6(e).  W  (L)  48-hr  prognosis  valid  OOZ  December  14,  1964. 
s 
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Isotach,  knots 
Error  field,  lO' 


Fig.  7(a).  Z(L)  12-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 
[(forecast  minus  observed)  over  areas  where  Wg(L)  observed  as  >  100  knots)]. 
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Fig.  7(b).  Z{L)  24-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  7(c).  Z(L)  36-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  8(a).  Z(L)  12-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 


Fig.  8(b).  Z(L)  24-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Isotach,  knots 
Error  field,  10^  ft 


S  / 

^^30*  — 

k>. 

Fig.  8(c).  Z(L)  36-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  8(d).  Z(L)  48-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  9(a),  12-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  9(b).  24-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  9(c).  36-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 


Isotach,  knota 


-3  ,  - 


- Error  field,  knot  10  ft 


Fig.  9(d).  48-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  10(a).  12-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 


Fig.  10(b).  24-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  10(c).  36-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  10(d).  48-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  11(a).  S  12-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  11(b).  Sa  24-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  11(c).  36-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  11(d).  Sa  48-hr  prognosis  (valid  OOZ  December  13,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  12(a).  Sg^  12-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig.  12(b).  24-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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Fig,  12(c).  Sa  36-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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- Error  field,  knot  10 
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Fig.  12(d).  Sa  48-hr  prognosis  (valid  OOZ  December  14,  1964)  error  field 

[(forecast  minus  observed)  over  areas  where  W  (L)  observed  as  >  100  knots)]. 

s 
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SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  feasibility  study  for  the  LMW  Predictor  Network  Technique  has  demon¬ 
strated  that  prognoses  of  12-  and  24-hr  LMW  wind  speed,  height,  shear  above  and 
shear  below  using  a  grid  oriented  according  to  the  LMW  wind  direction  are  reason¬ 
able  and  much  superior  to  persistence.  The  predictor  fields  comprising  the  initial- 
state  and  previous  12-hr  changes  of  LMW  wind  speed,  height,  shear  above  and  shear 
below,  and  the  initial-state  wind  speed  at  constant  pressure  surfaces  all  provide 
useful  prediction  information.  The  categorization  procedure  used  was,  in  general, 
successful.  The  separation  into  polar  and  subtropical  systems  and  the  classification, 
within  each  system,  of  grid  points  on  or  near  the  jet  core,  to  the  left  of  the  core,  and 
to  the  right  of  the  core,  proved  worthwhile.  The  attempt,  however,  to  further  categorize 
grid  points  along  the  core  by  their  location  relative  to  wind  speed  maxima  was  judged 
to  be  unnecessary.  This  result  is  perhaps  surprising  but  it  must  be  remembered 
that  the  horizontal  variations  of  the  predictands  along  the  jet  core  are  usually  smaller 
than  the  changes  experienced  to  the  left  and  right  of  the  core. 

While  the  feasibility  study  yielded  promising  results,  additional  effort  would 
be  required  to  develop  an  operationally-useful  prediction  technique.  An  expanded 
effort  should  be  characterized  by  the  following: 

(a)  larger  data  sample;  a  sample  containing  at  least  twice  as  many 
observation  times  as  that  used  in  the  feasibility  study  should  be  used  to 
insure  the  derivation  of  stable  relationships  for  the  subtropical  jet  stream, 

(b)  revision  of  grid-point  categorization  procedure;  the  three  categories 
along  the  jet  core  should  be  combined  into  a  single  category,  and  a  new 
category  containing  grid  points  further  downstream  from  the  initial-state 
position  of  the  jet-stream  wind-speed  maximum  should  be  added, 

(c)  use  of  prognostic  data;  operational  prognoses  of  constant-pressure- 
surface  height  and  wind  speed  should  be  included  in  further  screening-regression 
experi  ments . 

Evaluation  of  the  prognoses  of  LMW  parameters  generated  by  using  the  diagnostic 
LMW  equations  (originally  developed  for  generating  initial-guess  fields  for  the  LMW 
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analyses)  with  numerical  prognoses  from  the  six-level  dynamic  model  in  operation 
at  Offutt  AFB,  indicated  that  the  12-  and  24-hr  prognoses  were,  in  general,  a  good 
representation  of  the  observed  fields.  The  36-  and  48-hr  prognoses  deteriorated, 
but  were  still  considered  at  least  fair  and,  occasionally,  good  over  some  areas  of 
the  Northern  Hemisphere. 

It  is  concluded  that  the  quality  of  the  LMW  prognoses  generated  with  the 
diagnostic  LMW  equations  applied  to  constant-pressure-surface  prognoses  is  largely 
a  function  of  the  quality  of  the  constant -pressure-surface  prognoses.  Any  improve¬ 
ment  in  the  numerical  prognoses  will  be  reflected  by  an  improvement  in  the  LMW 
prognoses. 
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APPENDIX 

EQUATIONS  for  POLAR  and 
SUBTROPICAL  JET-STREAM  CATEGORIES 


Polar  Jet-stream  Categories 


1.  Category  lA 


W  (L)  =  17.095  +  0.4877[W  (2)(3,3)]  +  0.4487(W  (L)(4,  2)]  -  3.2824[AS  (7,  1)] 

S  S  S  D 

AZ(L)  =  1546.5  -  0.4406[Z(L)(5,3)]  -  0.2616[AZ(L)(5,  3)] 

1. 2 

^  1.1948  -  0.9698  [S  (5,3)]  +  0.0214[W  (L)(4,  2)] 

O  D  S 

AS  =  -  1.0189  -  0.8770[S  (5,3)]  -  0.0216[W  (L)(2,  2)] 

3.  3i  S 

W  (L)  =  18.290  4-  0.6189[W  (L)(3,3)]  +  0.3682[W  (2)(3,1)]  -  4.2949[S  (3,3)] 

s  s  so 

AZ(L)  =  1907.7  -  0.9647[Z(L)(5,3)]  +  0.3170[Z(L)(1,  3)]  +  6.4861[W  (1)(5,3)] 

S 

AS^24  =  1.4101  -  1.0273[S^(5,3)]  +  0.0266[W^(2)(3,  1)] 

AS  =  -  1.4198  -  0.9419[S  (5,3)]  -  0.0230[W  (2)(3,1)] 

8i  24  3,  S 


Category  IB 


AZ(L) 


12 


-  4.8861  +  0.3005[W  (L)(3,3)]  +  0.4266[W  (5)(5,3)]  +  0.3240[W  (L)(6,  2)] 

s  s  s 

+  0.2960[W  (1)(5,3)] 
s 

2059.5  -  0.9412[Z(L)(5,3)]  +  0.2901[Z(L)(3, 3)]  +  4.0053[W  (2)(3,1)] 

s 

-  4.7071[W  (5)(5,3)]  +  3.0909[W  (L)(5,  5)] 

s  s 


AS 


bl2 

iZ(L), 


'24 


AS 


b24 


=  1.1185  -  1.0498[S  (5,3)]  +  0.0259[W  (L)(6,  2)] 

D  S 

=  -  0.7059  -  0.9096[S  (5,3)]  -  0.0170[W  (3)(3,3)]  -  0.0194[W  (5)(5,3)] 

3.  S  S 

=  -  5.5470  +  0.3207[W  (L)(3,3)]  +  0.4771[W  (L)(6,  2)]  +  0.3189[W  (2)(3,5)] 

s  s  s 

=  2362.2  -  0.9523[Z(L)(5,3)]  +  0.2539[Z(L)(3,  1)]  -  0.1629[AZ(L)(7,  3)] 

=  0.9473  -  l.0868[S  (5,3)]  +  0.0177[W  (L)(6,  2)]  +  0.0177[W  (1)(3,  1)] 

Os  s 


AS 


a  24 


0.6210  -  1.0106[S  (5,3)]  -  0.0284[W  (2)(3,3)]  -  0.2396[S  (7, 1)] 
3.  S  0 


+  0.0235[W„(1)(3,3)]  -  0.0006[Z(L)(3,  1)] 
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3. 


Category  1C 

W  (L)  =  10.612  +  0.4848[W  (L)(5,3)]  +  0.2680[W  (L)(3,3)]  +  3.2712(8  (5,3)] 

S  12  S  S3, 

+  0.1915[W  (L)(4,  2)1 
s 

AZ(L)  =  1142.2  -  0.8010[Z(L)(5,3)]  +  0.2271[Z(L)(3, 3)]  +  0.1788[Z(L)(7, 3)] 

JL  ^ 

+  2.6849[W  (L)(5,l)] 
s 

AS.  _  =  0.6738  -  0.8300(8  (5,3)1  +  0.0203 [W  (L)(4,  2)] 

D  12  D  S 

AS  -  -  1.3579  -  0.8892(8  (5,3)]  -  0.0182[W  (3)(3,3)] 

3,  12  3,  S 

W  (L)^.  =  0.6702  +  0.3772[W  (1)(5,3)]  +  0.1464[W  (3)(1,3)]  +  0.2183[W  (L)(6,4)] 

s  24  s  s  s 

+  0.2100[W  (3)(3,3)]  +  0.2625[W  (5)(3,5)]  +  0.1562[W  (2)(3,1)] 
s  s  s 

AZ(L)  =  1576.1  -  0.9675[Z(L)(5,3)]  +  0.2901[Z(L)(7, 3)]  +  0.2366[Z(L)(3, 3)] 

AS.  =  0.8773  -  1.0092(8,  (5,3)]  +  0.0158[W  (L)(6,4)]  +  0.0115[W  (3)(1,3)] 

D  24  D  S  S 

AS^24  =  "  1-2374  -  0.9389(8^(5,3)]  -  0.0187[W^(L)(2, 4)]  +  0.0096[AW^(L)(3, 5)] 

4.  Category  ID 

W  (L)  =  2.2198  +  0.4372[W  (2)(5,3)]  +  0.3545[W  (L)(4,  2)]  +  0.2724(AW  (L)(3,3)] 

s  12  s  s  s 

+  0.3065[W  (L)(6,4)] 
s 

AZ(L),„  =  1761.2  -  0.71C8[Z(L)(5,3)]  +  0.2557[Z(L)(3, 3)]  -  4.0631[W  (3)(7,3)] 

12  S 

„  =  0.8290  -  1.0071(8,(5,3)]  +  0.0133[W  (L)(4,  2)]  +  0.0130[AW  (L)(3,3)] 

D 12  b  s  s 

+  0.0160[W  (2)(5,3)] 
s 

AS  =  -  0.6257  -  0.9887(8  (5,3)]  -  0.0120[W  (L)(4,  2)]  -  0.0154[AW  (L)(3,3)] 

3.  12  3,  S  S 

-  0.0216[W  (3)(5,3)] 
s 

W  (L)  =  21.646  +  0.3392[W  (2)(3,  3)]  +  0.1752[W  (L)(5,  1)]  +  0.2602[AWg(L)(3 , 3)] 

S  24  S 

+  0.3257[Wg(2)(5,3)] 

AZ(L)  =  1897.2  -  0.6691[Z(L)(5,3)]  +  8.1958[W  (1)(3,3)]  -  0.2088 [AZ(L)(5, 3)] 

24  S 
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AS, 


b24 


AS 


a  24 


1.7910  -  1.0590[S  (5,3)]  -  0.0186[AW  (L)(7,5)]  +  0.0135[W  {L)(4,  2)] 
D  S  S 

1.5516  -  0.8312[S^(5,3)]  -  0.0180[W^(2)(3, 3)] 


5.  Category  IE 

S'  ^12 


W  (L).„  =  3.9027  +  0.3656[W  (L)(5,3)]  +  0.2151[W  (L)(4,  2)]  +  0.2195[W  {L){7,5)] 


AZ(L) 


12 


AS, 


bl2 


^"al2 

^<"^^24 


AZ(L) 


24 


AS 


b24 


AS 


a  24 


+  0.1894[W  (3)(3,3)] 
s 

743.78  -  0.5463[Z(L)(5,3)]  +  0.3422[Z(L)(3, 3)]  -  0.1827 [AZ(L)(5,  3)] 
1.3318  -  0.8551[S  (5,3)1  +  0.0108[W  (L)(3,3)] 

D  S 

-  0.5632  -  0.7209[S  (5,3)]  +  0.1977[S  (7,5)]  -  0.0121[W  (3)(3,3)] 

3.  3.  S 

1.4805  +  0.3101[W  (L)(5,3)]  +  0.2819[W  (L)(2,4)]  +  0.3125[W  (L)(6,  2)] 
s  s  s 

-  3.3884[S  (7,5)] 

3. 

1521.4  -  0.8865[Z(L)(5,3)]  +  0.4084[Z(L)(3, 3)]  +  5.3421[W  (1)(7,3)] 

s 

1.6615  -  0.7991[S  (5,3)]  -  0.1966[S  (7,5)] 

D  3. 

-  0.5484  -  0.7733[S  (5,3)]  +  0.2505[S  (7,5)]  -  0.0137[W  (L)(l,  1)] 

3.  3.  S 
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Subtropical  Jet-stream  Categories 


6.  Category  2A 


W  (L) 
s 

AZ(L) 


12 

12 


AS 


bl2 


AS 


a  12 


AZ(L) 


24 


AS 


b24 


AS 


a  24 


1.5142  +  0.5737[W  (L)(5,3)]  +  0.3287[W  (3)(5,3)] 
s  s 

3805.6  -  0.8121[Z(L)(5,3)]  -  12.520[W  (5)(3,3)]  +  6.9727[W  (1)(1,3)] 

s  s 

4.4566[W  (L)(7,3)]  +  0.2311[AZ(L)(1, 3)]  +  3.7162[AW  (L)(7,  1)] 
s  s 

2.8210  -  0.8550[S^(5,3)] 

0.1014  -  0.8701[S  (5,3)]  -  0.0330[W  (3)(5,3)]  -  0.1896[AS  (5,3)] 

3.  S  D 

29.650  +  0.7894[W  (L)(6,  2)] 
s 

3887.3  -  0.9214[Z(L)(5,3)]  -  6.3240[W  (3)(3,3)]  -  50.623[S  (7,3)] 

S  3 

3.5424  -  0.8072[S  (5,3)]  -  0.0332[AW  (L)(5,3)] 

D  S 

2.8107  -  0.8164[S  (5,3)] 

3 


7 .  Category  2B 

W  (L)  =  -  0.8421  +  0.5364[W  (L)(5,3)]  +  0.4871[W  (L)(4,  2)] 

S  S  S 

AZ(L)  =  3093.4  -  0.6087[Z(L)(5,  3)]  -  7.8424[W  (L)(2,  2)] 

s 


AS 


bl2 


=  -  0.9554  -  0.3656[AS  (5,3)]  +  0.0301[AW  (L)(l,3)]  -  0.0216[AW  (L)(7,5)] 

D  S  S 


0.3157[S  (5,3)]  +  0.0182[W  (L)(4,  2)] 
o  s 


W  (L) 


24 


1.2432  -  0.5327[S  (5,3)]  -  0.0308[W  (2)(5, 3)] 

3  S 

20.087  +  0.3186[W  (L)(4,  2)]  +  0.3224[W  (L)(5,3)]  +  4.7280[AS  (3,  1)] 
s  s  o 


-  3.9829[S^(3,3)] 


AZ(L) 


24 


=  2785.5  -  0.7130[Z(L)(5,3)]  -  73.241[AS^(7,  5)]  -  68.015[AS^(3,  1)] 


-  60.323[S  (7,3)]  +  48.199[S  (5,3)] 

3  3 


AS 


b24 


0.7058  -  0.7571[S^(5,3)]  -  0,4444[S^(3,  1)] 


AS 


a  24 


0.6807  -  0.5762[S  (5,3)]  -  0.0263[W  (L)(4,  2)] 

3  S 
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8.  Category  2C 


W  (L)  =  12.144  +  0.4928[W  (L)(5,3)]  +  0.4321[W  (L)(5,  1)] 

S  12  S  S 

AZ(L)  =  2177.2  -  0.2932[AZ(L)(5,3)]  -  11.112[W  (5)(3,3)]  -  0.4261[Z(L)(5, 3)] 

s 

-  3.4979 [AW  (L)(3,5)] 

s 

=  1.8988  -  0.6107[S  (5,3)]  +  0.4323[AS  (7,3)] 

D  12  D  3, 

AS  =  -  1.4245  -  0.7990[S  (5,3)]  +  0.4787[S  (5,1)] 

3  12  3  3 

W  (L)  =  -  3.7601  +  0.3850[W  (2)(5,3)]  +  0.7078[W  (L)(5, 1)]  +  0.2830[W  (L)(3,5)] 

S  24  S  S  S 

-  6.4601[S^(5,1)] 

AZ(L)  =  2240.5  -  0.5834[Z(L)(5, 3)] 

ASj^24  =  "  0-2932  -  1.0465[S^(5,  3)]  +  0.0231[W^(2)(5,  3)]  -  0.3866[AS^(1,  3)] 

+  0.0017[AZ(L)(3,3)]  +  0.5484[AS  (1,1)]  +  0.0315[W  (3)(3,1)] 

3  S 

AS  =  -  2.7770  -  0.8112[S  (5,3)] 

3  24  3 


9.  Category  2D 


W  (L)  =  12.141  +  0.4560[W  (L)(5,3)]  +  0.3706[W  (2)(3,3)]  +  0.2835[W  (5)(5,3)] 

S  12  s  s  s 


AZ(L) 


12 


12 


506.6  -  0.5185[Z(L)(5,3)]  +  0.3897[Z(L)(3, 3)] 
1.1400  -  0.6452[S^(5,3)]  +  0.3062[S^(3, 3)] 


AS 


a  12 


0.4872  -  0.8305[S  (5,3)]  -  0.0252[W  (L)(3,3)] 

3  S 


3.8416  +  0.6138[W  (L)(4,  2)]  +  0.3972[W  (2)(3,3)] 
s  s 


AZ(L) 


24 


3033.6  -  0.7333[Z(L)(5,3)]  -  7.000[W  (5)(1,3)] 

s 


AS, 


b24 


0.2924  -  0.9657[S  (5,3)]  +  0.0281[W  (L)(4,  2)] 
0  s 


AS 


a  24 


0.6528  -  0.9875[S  (5,3)]  -  0.0265[W  (L)(4,  2)] 

3  S 
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10.  Category  2E 

W  (L)  =  15.975  +  0.4418[W  (3)(5,3)]  +  0.4946[W  (1)(7,3)]  +  0.3421[W  (5)(3,3)] 

S  S  -  s  s 


AZ(L) 


12 


3582.2  -  1.0579[Z(L)(5,3)]  -  6.5998[W  (3)(5,3)]  +  0.2748 [Z(L)(3 , 3)] 

s 


7.5708[W  (5)(5,3)]  +  4.5680[W  (L)(3,3)] 
s  s 


no 

b  12 

=  1.9054 

=  -  1.2837 

=  16.673 

AZ(L)24 

=  3724.2 

'^b24 

=  2.0674 

^a24 

=  -  1.0445 
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